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The  prioary  objective  of  this  study  »as  to  provide 
a co.otribution  to  the  phenonecon  "Suppressioa"  as  an 
aspect  vithin  the  military  envixonment. 

Analytical  models  explaining  these  aspects  were 
developed  in  order  to  identify  the  influences  to 
suppression.  Techniques  are  examined  for  including  the 
suppressive  effects  of  weapon  systems  in  Lanchester 
type  combat  models,  which  may  be  useful  in  wargame 
evaluations  of  military  judgements,  and  in  force  level 
planning . The'  study  also  provides  techniques  to 
analyze  and  fit  experimental  data  to  the  analytical 
models. 

The  data  to  verify  the  models  were  obtained  from 
related  experiments  performed  by  Combat  Development 
Experimentation  Command  (CDEC) , Fort  Ord,  California. 

The  result  for  the  modelling  approach  to 
suppression  indicates  source  dependences  on 
quantitative  as  well  as  on  qualitative  features. 

The  functions  are  left  quite  general,  although 
some  functional  forms  are  derived  and  discussed. 
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I*  gORMOLATION  OF  PROEIEM 


a.  lETBOOOCIION 

In  recent  years  considerations  of  Human  Factors  aspects 
in  the  military  has  gained  more  and  more  importance.  The 
study  of  psychological,  physiological,  and  environmental 
conditions  and  their  influences  on  the  performance  of  men  in 
“Man-ilachine-systems**  supports  the  development  of  new 
doctrines,  design  of  weapon  systems  as  well  as  training 
programs  for  troops. 

One  phenomenon  in  military  man*-machine-systems  is 
suppression.  Modelling  this  phenomenon  has  only  recently 
been  given  much  attention. 

This  is  partly  because  modellers  did  not  understand 
exactly  neither  the  causes  of  suppression  nor  how  it  affects 
the  course  of  combats. 

There  is  an  intuitive  feeling  that  when  a soldier  or  a 
combat  unit  is  being  fired  upon,  it  will  be  less  effective 
than  when  it  is  not  receiving  fire.  This  is  generally 
referred  to  as  suppression,  but  it  can  include  much  more.  In 
the  broadest  sense,  suppression  can  effect  individuals, 
units,  or  weapon  systems  in  different  types  cf  combat. 

This  paper  uill  limit  itself  to  individuals  or  small 
infantry  units. 


7 


In  order  to  be  able  to  build  up  functional  relationships 
based  on  the  suppression  idea^  between  different  categories 
of  people  and  tioe,  this  paper  will  first  exanine  different 
rates  as  input  to  analytical  Lanchester-type-oodels.  The 
rates  will  be  based  on  stochastic  aspects.  The  sodelllng  is 
done  under  different  viewpoints;  their  results  will  be 
ccspared  and  discussed. 

In  the  second  part  of  the  thesis#  the  parameters  cf  the 
■odels  will  be  tested  in  a regression  analysis  against  real 
world  data  which  were  placed  at  the  writers'  disposal  by 
CDEC, 


The  goal  of  this  paper  is  not  tc  present  a final 
fraoeworJc  cf  suppression;  however  iz  may  contribute  to 
clarification  of  aspects  of  suppression#  and  help  to  embed 
it  as  a component  in  future  large  size  models. 


B.  BACKGBOQNO 


The  strains  imposed  on  individuals  in  our  society  are 
constantly  increasing.  Modern  technologies  and  constant 
efforts  for  improvenents  of  standards  of  living  lead  to 
growing  difficulties  in  adjustments  or  sometimes  to  complete 
failures  to  adjust.  These  facts  create  stress  and  we  may 
observe  that  the  degree  of  stress  increases  with  the 
difficulty  of  tbe  ad justaent-problea. 

The  tern  stress  will  be  used  as  a substitute  for  what 
light  be  called  otherwise  as  anxiety,  conflict#  eiotional 
distress#  extrene  environoental  conditions#  ego-threat, 
frustration#  threat  to  security,  tension  or  probably 
arousal.  Stress  can  be  thought  of  as  the  result  of  aliost 
any  envircniental  interference  £ Appley#p.2  ]. 


8 


The  stress^generatlng  . features  of  the  civilian 
environaent  are  great,  hut  the  environnenr  created  by  modern 
warfare  possesses  additional  features  which  result  in  an 
increase  of  stress.  The  coobat  environaent  created  by  the 
weapon  power  of  the  enemy  causes  a constant  threat  to  life. 
The  soldier  has  to  operate  under  this  threat  and  naturally 
he  will  respond  with  constantly  recurring  fear.  This  may 
break  down  the  soldier's  psychological  and  physiological 
resistance. 

. Fear  and  anxiety  in  battle  is  common,  being  experienced 
by  between  80  and  90  % of  combatants.  Pains  in  the  stomach, 
fatigue,  dizziness,  perspiration,  and  enhanced  heart-beat 
are  some  vegetative  correlates  to  fear  and  an.\:iety.  Of 
course  the  moment  when  the  individual  soldier  reaches  his 
breaking  point  raries  and  depends  cn  individual 
psychological  and  physiological  resistance  and  the  severity 
of  the  battle. 

One  interesting  observation  from  Agiell  was  that 
auditory  sensations  corvey  the  stress  of  battle  most 
strongly  and  most  directly.  The  psychological  effect  of 
weapons  is  directly  related  tc  their  sound  level  and  the 
frequency  with  which  the  sound  occurs  [ Agrell,p. 215 ]. 

Each  enemy  grenade  causes  the  soldier  to  react 
constantly  with  fear.  Stress  and  fear  can  have  a 
significant  sensory  operating  characteristic,  e.  g.  the 
detection  threshold  and/or  sensitivity  may  decrease  as  a 
result  of  stress  (Weltoann,  6,  Christianson,  3. A.  and 
Egstroa,  G. H.,  p. 423*430] • 
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C.  tULlTAai  SOPPfi-SSSXCli  IN  COBBAS-ESVXBOilHEBI 


When  weapons  are  used  in  combat,  there  are  two  types  of 
effects  that  they  have.  The  first  type  of  effect  is  physical 
damage  or  injury  to  the  target  and  the  second  type  of  effect 
is  psychological. 

This  second  effect  of  a weapon  has  led  to  the  term  or 
concept  “Combat  Suppression”.  There  exist  a general  belief 
that  fire  suppression  is  important,  but  the  importance  of 
suppression  effect  on  combat  outcoimes  as  compared  to  the 
effect  of  other  areas  such  as  firepower,  mobility, 
intelligence,  command/control  has  not  been  quantified 
adequately. 

1 • Definitions 

Suppression  can  be  generally  defined  as 
the  temporary  degradation  in  the  quality 
of  performance  of  an  individual  soldier 
or  unit  by  an  internal  or  external 
stimulus. 

CDEC,  1977 

A more  useful  operational  definition  in  terns  of 
performance  capability  changes  is  provided  by  the  "Repcrt  of 
the  Army  scientific  Adviscry  panel  Ad  hoc  Group  on  fire 
Suppression”. 


Fire  suppression  is  a process  which 
causes  temporary  changes  in  performance 
capabilities  of  the  suppresses 
{suppressed  system)  from  those  expected 
when  functioning  in  an  environment  he 
Icnows  to  be  passive.  Thes^  changes  are 
caused  by  signals  from  delivered  fire  or 
the  threat  of  delivered  fire,  and  they 
result  from  behaviors  tnat  are  iotondea 
to  lessen  risk  to  the  suppressee. 

Ad  Hoc  Group,  1976 

This  definition  emphasizes  that  suppcesslon  is  not  a 
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single  effect  which  can  be  oeasured  totally  on  a single 
quantitive  scale.  Suppression  effects  ace  Dultidimensional 
and  the  "amount  of  suppressicn"  varies  among  these 
dimensions  (e.g.  fire  impact  points,  soldier's 
charactecists,  and  reaction  to  the  fire,  his  combat 
experience  etc.}* 

Suppressive  fire  in  a combat  environment  can  supp'ress  a 
number  of  combat  activities;  for  example:  firing,  search  fo/ 
and  observation  of  targets,  movements  of  units  or  ccmmar.'4 
and  control. 

2*  Structure  0|  Fire  Su ppression  Process 

lire  suppression  is  a complicated  process  involving 
many  physical,  environmental,  physiological,  behavioral,  and 
operational  variables.  The  important  point  to  emphasize  is 
that  the  behavior  involved  is  in  response  tc  stimuli  that 
originate  both  externally  (combat  environment}  and 
internally  (personal  background,  training  and  experience)  to 
the  soldier  suppressee.  These  aspects,  however,  are  not 
iucluded  in  this  paper. 

The  intensity  and  duration  of  suppression  can  not  be 
predicted  from  a knowledge  of  the  combat  environment  alone. 
It  requires  an  analysis  of  the  underlying  motivational  and 


cultural  factors 

environment. 

and 

of 

the  context  of 

the  combat 

The  Fig  1 

C hd 

Hoc 

Group,  p.  36  ] shows 

a schematic 

desccipticn  of  a process  when  suppressive  fire  is  delivered 
and  its  affect  on  the  combat  result. 


the  process  is  discribed  b^  individual  functicns* 


a.  Signal  Process 

Ihe  first  process  in  fire  suppression  is  the  genecaticn  of 
signals  provided  by  suppressing  weapons. 

Inputs  to  this  process  are  the 

characteristics  of  the  weapon  systems  (caliber « 
amount  of  porpellant,  warheadtypes  etc)  and 

the  environmental  characteristics  (trajectory, 

platform,  arrival  points  etc). 

Characteristics  of  the  weapon  systems  can  vary  in  order  to 
increase  or  decrease  suppressive  effects. 

Some  parameters  that  are  considered  to  be  important  to 
suppression  signals  are: 

huzzle  velocity  (an  increase  in  muzzle  velocity  is 
associated  with  an  increase  in  signal  variables)  ; 

Caliber  (as  caliber  increases,  the  firing  signals 
and  projectile  signals  increase  along  with 
lethality)  ; 

Projectile  weight  ( penetration  depends  on  weight 
and  velocity  at  impact  and  increases  shoch  coupling 
to  ground)  ; 

Harhead  charge  weight  ( the  explosive  charge  weight 
determines  the  energy  in  the  pressure  pulse) ; 

Additional  parameters  like  fire  frequency  and 
proximity  of  shots  could  also  be  mentioned. 

Snvironmental  characteristics  are  also  variables  but  they 
can  not  be  determined  completely. 

Environment  has  an  influence  on  the  signal 
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generation  -and  transraission-process.  for  exaapXe^  auditory 
signals  that  result  £roa  the  impact  o£  projectiles  depend 
heavily  on  the  nature  o£  the  object,  or  material  impacted* 

A soft  yielding  material  such  as  dusty  ground  gir 
sand  receiving  the  impact  of  a projectile  will  produce  a 
different  pulse  and  sound  than  will  hard  unyielding  ground 
under  the  sane  impact. 

Sound  signals  can  be  attenuated  by  the  shadowing 
effect  of  large  obstacles  or  may  be  increased  by  e'rho  or 
reverberation.  Visual  signals  are  strongly  modified  by  the 
condition  of  lighting.  Haze,  fog,  rain,  and  snow  act 
similarly  to  smoke  and  dust.  The  visual  field  is  also 
reduced  and  interrupted  by  terrain  and  other  obstacles. 

The  envlrcument  modifies  the  produced  stimuli 
when  they  are  transmitted  to  the  location  of  the  suppressee. 
Outputs  of  this  process  are  the  attenuated  sensory  signals 
that  become  imputs  to  the  human  sensory  receptors. 

b.  Human  Process 

Hany  of  the  determinants  of  the  soldier's  performance  on 
the  battle  field  are  unkown  or  at  least  uncertain  - thought 
of  as  influenced  by  chance  factors.  This  emphasizes  the 
difficulties  of  predicting  human  behavior  in  a combat 
environment.  Tae  human  process  (sensory  and  perception) 
converts  the  received  signals  into  a perception  of  the  risk. 

Eattie field  stimuli  effecting  the  individual  are 
detected  and  converted  into  sensory  data  by  a process  such 
as  vision  and  audition,  so  the  sensory  process  suggests  that 
the  weapon  systems  stimuli  relevant  to  suppression  are  the 
-loudness,  and 
-visual  impact. 

There  exist  moderating  factors  that  influonoa 
the  operating  char act eristcs  of  the  sensory  process  and  that 
determines  which  stimuli  ace  effective, 
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Sensory  nodifiers  (i.e«  earplugs# 

eight-vis loD»devices)  serve  to  change  .users*  sensitivity 
range.  A oajor  effect  of  these  devices  is  to  change  the 
salience  of  stiouli. 

High  concentration  on  an  activity  or  a high 
level  of  effort  on  an  activity  {e.g.  oisslle-gunner  is 
tracking  a target  or  reloading  bis  system)  nay  increa-se  the 
absolute  thresholds 

Ihe  posture  of  a soldier  (standing  or  sitting) 
and  the  sequence  of  posture  (observing,  ducking,  observing) 
influences  the  sensory  capabilities  (e.g.  observing  for  10 
seconds  continuously  is  not  equivalent  to  observing  5 

seconds  , ducking  10  seconds  and  observing  6 seconds) . 

Ihe  perception  process  integrates  sensory  and 
other  information  into  a perception  of  the  risk.  Risk  refers 
to  the  uncertainty  of  damage,  injury,  or  loss.  It 

char'vcterizes  decision  situations  in  which  the  consequences 
of  choosing  an  action  are  uncertain. 

If  there  is  no  uncertainty  in  the  possible 

wOtcoae,  there  is  no  risk. 

Esreeived  risk  is  a function  of  uncertainty  and 
the  subjective  value  the  individual  associates  with  each 
outcome;. 

Ferceived  risk  represents  the  output  of  the 

combined  ’senso^-y  and  perception  process.  It  depends  en  the 
individual's  experiences  and  training  in  assessing  risk  from 
sensor?  infurmation,  Also  cover  provided  by  the  environment 
and  the  individual's  posture  iB?.y  influence  risk  perception. 

c,  ^Beact^ion  Process 

aiveu  the  input  perception  of  risk,  this  process  causes 
physical  and  mental  reactions,  which  depend  on  the 
-current  mission 
-task 
^activity 


*coabat  training  doctrine  and  experience 

-group  dynaelc^i  and 

-the  guallty  o£  leadership. 

It  is  conjectured  that  tvo  individuals  vho  perceive  the  sane 
degree  of  high  risk  , but  who  have  different  amounts  of 
ccnbat  engageaent  experience,  night  be  likely  to  react 
differently  tc  the  risks. ' 

The  soldier's  reaction  is  also  influenced  by  his 
curreU;t  state.  A soldier  who  has  recently  ducked  nay  be  noce 
likely  to  duck  than  one  who  has  not  , given  the  sane 
delivered  fire. 

Prior  reaction  or  sequence  of  reactions  may  be  a 
good  predictor  of  the  coning  reaction. 


d.  Performance  Effects  Process 

Given  the  reactions  of  the  hunan  behavior  process,  .it  is 
conjectured  that  these  directly  affect  the  performance  of 
certain  activities  of  the  suppressee  in  a calculable  way. 

If  for  exasple  the  suppressee  takes  cover,  he 
nay  fire  less  often  and  less  accurately  and  also  night  be 
less  vulnerable.  The  nagnitude  and  duration  of-  these  changes 
in  performance  are  dependent  on  the  characteristics  of  the 
systen  employed  by  the  suppressee  and  the  target  of  his 
activity. 

So  the  nature  and  duration  of  change  in 
perfornance  capabilities  is  deternined  by  the  performance 
effects  process. 


3.  .Su£££a£S42B  la  llaM 

A fire  suppression  research  prograa  requires 
significant  experinentation  on  behavioral  attitudes  and 
reactions  tc  risk.  This  necessity  causes  treaendous 
difficulties  in  trying  to  induce  actual  behavior  in  soldiers 
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in  field  experinents.  Forner  studies  sbcued,  that  the 
soldiers  felt  true  psychological  stress  only  in  situations 
in  which  they  believed  that  they  were  in  real  danger.  Such 
situations  are  difficult  to  contrive  and  to  control.  Social 
and  ethical  liaits  and  legal  regulations  preclude  the 
introduction  of  actual  physical  risk,  soldiers  niist  be 
taught  the  “rules  and  risk”  defined  in  that  context'.  The 
success  of  playing  the  role,  being  an  individual 
participating  in  a conbat  engagement,  depends  on  the 
soldiers'  motivation  and  willingness. 

Because  of  these  reasons  and  the  sultidimensional 
shape  of  the  fire  suppression  process  as  mentioned  earlier, 
suppression  in  field  experimentations  may  be  restricted  only 
to  some  variables  involved  in  this  process. 

D.  APPaOPBJATB  OBJECTIVES 

The  overall  objective  of  a fire  suppression  research 
should  be  to  relate  changes  in  performance  capabilities 
caused  by  fire  suppression  [Ad  l^oc  Group,  p.  110]. 
Besponding  more  directly  the  following  objectives  may  be 
determined: 

Indicating  the  effects  of  suppression  on  combat 
results,  i.  e.  to  develop  rates  cf  suppression. 
These  values  nay  be  compared  to  other  effected  areas 
and  probably  employed  in  computer  sixulations.  These 
numbers  represent  two  kinds  of  variables:  Weapon 
system  variables  and  human  suppression  performance, 
given  operational  and  environmental  conditions. 

Cetermining  characteristics  cf  suppressive  fire 
systems,  characteristics  which  should  be  assigned  to 
such  a weapon  system.  Besults  are  developed 
experimentally  . Chapter  ill  B.  will  display  some 
evaluated  parameters  and  constants  for  the  developed 
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aicdel,  which  nay  also  repraseat  the  supprassive 
characteristics  o£  the  used  weapon  systens* 

Baducing  suppressive  effects.  Hays  to  reduce  the 
effect  of  suppressive  fire  nay  also  be  considered  as 
an  appropriate  objective  of  suppression  research. 
Special  training  or  eguipnent  can  be  assigned  to  the 
soldiers  or  new  tactics  can  be  developed.  This 
objective  is  beyond  of  the  research  of  this  paper. 

E.  POSSIBLE  ALIEBNATIVES 

In  order  to  get  infornation  about  the  fire  suppression 
process,  previous  investigations  were  based  on  interviews 
and  guestionaires,  because  valuable  infornation  of  the  fire 
suppression  process  is  stored  in  the  ninds  of  conbat 
veterans.  Studies  on  veterans  of  the  Vietnam  conflict  and 
the  wars  cf  the  Near  East  would  be  especially  useful,  since 
newer  weapons  were  employed  and  the  combats  were  shorter  and 
more  intensive. 

These  studies  may  provide  a good  insight  to  the 
suppression  process  and/or  may  also  deliver  valuable  inputs 
for  the  mcdelling  approach. 


II.  CONSTROCTIOH  OF  ^ ANALYTICAL  MODEL 


S; 


A.  ASSUMPTIONS  ANC  SIMPLIFICATIONS 

Based  on  the  foregoing  discussion^  the  following  model 
is  a detailed  model  [ Taylor , 1978, p. 12  ] which  starts  out  by 
considering  the  behavior  aspect  of  a human  being  under  the 
influence  of  artillery  fire  power.  It  is  assumed  that  such 
simple  models,  which  represent  a small  part  of  a total 
scenario,  can  be  used  profitably  to  investigate  system 
dynamics  of  more  complex  models.  The  value  of  the  model 
derives  from  the  fact  that  it  forms  intuitively  plausible 
and  transparent  subsets  in  a large  compcsition  of  other 
subsets  which  determine  the  basic  structure  of  the  complex 
operational  model.  In  other  words,  the  whole  is  described 
in  terms  of  the  sum  of  its  parts. 

The  basic  concern  of  the  analytical  model  developed  here 
will  be  to  model  the  behavior  of  an  individual  experiencing 
artillery  fire,  considered  as  a function  of  tine,  where  that 
behavior  depends  upon  ammunition  types  and  location  of 
detonating  rounds. 

On  the  basis  of  particular  assumptions  and 
simplifications  it  will  be  possible  to  apply  the  results 
obtained  to  a group  of  people  (a  force)  on  a battelfield. 
The  result  of  these  considerations  will  provide  a 
relationship  between  time  and  the  actual  number  of  people 
affected  by  the  fire  power  of  the  artillery  . This  last  step 
of  the  model  is  carried  out  by  using  Lanchestei  type 
eguations,  so  called  after  the  pioneering  work  of  F.  N. 
Lanchester.  Finally,  the  models  enable  one  to  estimate  the 
total  firepower  of  the  force  at  any  point  in  time;  depletion 
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o£  total  firepower 
suppcessiott. 


is  . caused  by  attrition  and  by 


E.  HATfiS  FOB  THE  MODEL 

* Rates  Of  Suppression 

The  basic  considerations  in  the  preceeding  paragraph 
support  the  assu  nipt  ions,  that  the  behavior  of  a suppresses 
can  be  expressed  by  a conditional  probability  of  suppression 
as  a function  of  oiss  distance  r and  aspect  angle  e.  This 
function  is  represented  by  the  fanily  of  surfaces  of  the 
foro: 

P(S/0,r)  = exp[-  •“  r^(l  - e cos  9)]  (2.1) 


where  P(S/6,r)  is  the  conditional  probability  of  suppression 
given  that  a parti cc'tr  round  iapacts  under  a certain  aspect 
angle  6 and  a certain  miss  distance  r away  from  the  foxhole. 

Ihe  line  along  the  angle  9=0®  is  identical  with  the 
line  of  sight.  It  is  the  nain  direction  of  observation. 

The  constant  K and  e are  paraaeters,  which  are 
deteroined  by  the  experiaent  itself  and  by  the  environoental 
conditions. 

They  can  be  influenced  by  factors  as  discussed  in 
chapter  I which  aay  be  recalled  here  briefly. 

-type  of  aoaunition 
-freguency  of  arrival  of  rounds 
-perceptual  daaage 

-total  tiae  spent  in  the  foxhole  (learning  process) 
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-noise  appearance  of  the  rounds 
-flash  light  intensity  of  the  rounds 
-perforoance 

-personal  factors  like  age,  personal  condition,  etc. 

-degree  of  stress 

-notivation 

The  conputational  evaluation  of  the  constant  K and  e will  be 
performed  in  chapter  III  B.  and  C.  In  particular  it  vill  be 
important  to  determine  K as  it  varies  vith  different  types 
of  amauniticn. 

The  mathematical  conditions  for  the  two  parameters  K and  e 
are: 

0 < e < 1 
K > 0 

If  9 is  held  fixed  ,0^0  060*  and  P(S/9',r)  varies 
between  O<P(S/0,r)<1  we  obtain  a family  of  functions,  which 
is  two-dimensional  and  shows  an  exponential  relationship 
between  P(S/6,r)  and  r. 

This  is  illustrated  in  the  following  figure,  where 
the  angle  6 is  held  fixed  at  C*and  180*  The  parameters  K and 
c are  assumed  as  being  1500  and  0.7  respectively. 
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for 

foxhole 

for  9,"0* 

Figure  2 - FUNCXICN 

1 Z, 

’?(S/9Q,r)  - oxp(- 

- e cos  9)] 

22- 
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f.  ? 
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If  CQ  the  other  hand  P(S/'0/i)  is  fixed, 
0<P (S/0,r)<1 , and  the  equation  is  solved  for  r,  ve  will 
obtain  an  "egg-shaped"  function  with  iso-levels  of 
probability  of  suppression  P(S/6,r).  The  foxhole  is  located 
in  the  middle  of  the  coordinate  system.  Along  the  range 
r is  a max  for  a certain  fixed  probability  cf  suppression  , 
while  at  6=180*  r is  a min  for  the  same  probability* 

The  function  takes  the  form: 


2 ^ -K  an  P(S/0.r) 
^ 1 - e cos  0 


(2.2) 


i 


The  following  graph  shows  the  family  of  functions  for  4 
representative  selective  probabilities  of  suppressions. 


P{S/0,r)=O.1‘i  i=1,2,3,4  % (2.3) 

The  parameters  K and  e are  again  assumed  to  be  1500  and  0.7 
respectively* 


c.- 
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Ihese  '‘egg-shaped»-functions  of  iso-pcobabilities  of 
suppressioo  sioulate  the  reaction  of  haaan  beings  looking 
along  the  nain  axis  6=0^ in  a very  simple  way.  The  functions 
take  into  account  the  visual  and  accoustical  perception 
resulting  from  any  given  detonation  of  a round,  where  noise 
and  light  are  the  major  stimuli. 

One  could  think  of  other  functions  wich  simulate  the 
behavior  of  an  antitank  gunner  exposed  to  artillery  rounds 
for  instance : 


1 2 

P(S/e,r)  « L*exp(-  ^ r (1  - e cos  6)] 


(2.4) 


Ihis  type  or  function  allows  the  probability  of  suppression 
E(S/0,r)  to  te  smaller  than  one  at  its  maxisai  value  and  has 
the  same  general  behavior  as  the  function  before.  &s  a 
third  modification: 


P(S/0,r)  » exp[-  ^M(l  ” e cos  9)] 


(2.5) 


these  functions  have  the  disadvantage  that  they  have  a 
discontinuity  at  r-O.  The  integration  which  is  necessary  in 
the  following  derivation  is  more  difficult  than  the  chosen 


However,  the  crosscut  sections  of  these  functions 
are  not  "egg-shaped"  iso- functions  but  rather  simple  conic 
sections  (ellipses)  where  the  foxhole  is  located  In  the 
center  of  one  focus  point, 

the  paper  will  continue  with  the  function  first 
described  in  (2.1)^  because  of  its  simplicity  and  variety 
cf  applicatlcn. 

In  order  to  derive  a rate  of  suppression  for  the 
model/  it  is  necessary  to  evaluate  now,  in  a second  step. 
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the  uncoQditicnal  pcobability  of  suppression  P(S)  for  any 
inconing  round,  no  natter  where  it  will  iopact  around  the 
foxhole,  this  will  be  performed  by  oatching  the  conditional 
probability  P(S/0,r}  with  the  area-hir-probakility  P(A), 

Consequently,  to  the  engagement  procedures  of  the 
artillery,  the  targets  which  shall  be  suppressed  by  the 
artillery  are  categorized  as  snail  personal  targets.'  This 
inplies  that  the  mean  point  of  impact  (MPI)  of  a given  sat 
of  rounds  lies  on  the  target,  which  means  also,  that  the 
density  has  its  max  value  at  that  point.  According  to  the 
Q.S,  Array  Field  Manual  FM  6-161-1  page  2-2  it  can  be 
ascertained  that  the  MPI-error  is  destributed  normally  with 
its  mean  at  the  aim-point.  Thus  it  will  be  hypothesised  that 
in  such  a case  the  distribution  of  incoming  rounds  is 
bivariate  nctmal  with  parameters 


.1  = 0 
.2  - 0 
**  *^2  “ ^ 


(2.6) 


In  order  to  simplify  the  model,  it  is  proposed  that  the 
dispersion  of  rounds  expressed  in  the  standard  deviation  is 
equal  in  ail  directions.  With 

P a 0 {2.7) 

20  correlation  is  assumed  between  horizontal  and  vertical 
deviation.  Ihe  normality  is  also  preserved  if  mote  than  one 
artillery  gun  is  shooting.  The  essential  change  which  has 
to  be  made  when  a whole  artillery  unit  will  deliver  the 
rounds  will  be  the  value  of  the  standard  deviation  cr 
In  addition,  a is  determined  by  factors  like; 

dispersion 

the  fact  that  the  location  of  the  target  is  only 
estimated  and  an  artillery  unit  delivers  rounds  in  a 
fire  are?..  . 


type  of  anounitlon 


conditions  of  the  weapon  systeas  (precision) 
ballistic  properties  of  the  rounds 


distance  to  target 


caliber 


wind  and  other  weather  conditions. 


It  follows  frcffl  the  bivariate  normal  density  for  artillery 
bitsj 


f(x,  ,x-)  5 — 5 

^ ^ 25r(l-p^)a^ 


(2.8) 


and  by  inserting  the  previous  mentioned  assumptions  and 
chnngiiig  to  polar  coordinates: 


2 2^2 
r - + X2 


(2.9) 


cos  6 = — 
r 


that  the  unconditional  probability  of  suppression  P(S)  can 
he  written  in  the  form: 


0*2ir  r*” 

P(S)  » / / P(S/0,r)*f (r,0)  derdr 

0»O  r“0 


(2.10) 


Becall  that  an  area  element  in  polar  coordinates  can  be 
expressed: 


r*d0*dr  ■ 


Ihe  conditional  probability  of  suppression  is 


P([S/8,r)  - exp(-  ^ r' (1  - oof  0)J 


(2.11) 


and  the  density  for  artillery  hits  f(r#6)  is 


f(r,0)  >•  expi-  j (~)^] 

2vc  2 ° 


inserting  both  into  the  expression  above  «e  can 
the  integration: 


0=2ff  r=«“  17-  1 1 2 

P(s)  = f / exp[-  - r (1  - G cos  6)  r exp[-  jih  ]• 

0=0  r=0  2ita 


After  a change  of  integration  variable; 

2 


z * 0 
r*»0 


d(-^)  a dz 


Z ^ ■ 09 

r w 


we  find 

P(S) 


0»2Tr  zao*  r / 9 1 \T  1 

/ / exp  I - (1  - e cos  0)  + ~ ) I . — 

0»O  z«0  /J  7^r 


2iro‘ 


0«2v 

PCS)  » / 


S»0  -2vo^  I ~ (1  - e cos  0)  + j 

[-(I  a -sees  n 


*exp 


0-2v 

PCS)  - / 


0-0  2vo^||-^cos  0 + --^ 


0 ' jz«0' 

dO 


d6 


ZUla 


■) 


PCS) 


0«2ff 


/ 


Auct  0«»O  /l  "f  + (“C)  cos  0 


>d0 


(2.12) 


perforn 

r2 

dC-y)  d0 
C2,13) 


C2.14) 


dz  d0 

C2.15) 


(2.16) 


by  teplitvcing 


b « 1 + - j . 
2o 


c ■ -e 


the  integral  reduces  to  a known  form  which  can  be  solved 


e^air 

Ks)  » 2 / ]3  + C*c 

4va 


cos  8 


2 V 2 

where  b > c 


,(S,  K Y^e ^ K 

2jto^  sio  ® 2iro^’  yX7  2 


(b-c)  tan 


^ It 


z?:? 


2 / 5 2 

2a  vb  -c 


Inserting  back  the  expressions  for  b and  c,  one  cbtains 


After  the  data  analysis  in  chapter  IXI  B.  and  C. 
where  the  paraeeters  K#  e and  cr  will  be  deteroinedf  it 
will  be  possible  to  evaluate  the  probability  of  supi^ession 


e(S}  accoiding  to  the  proceeding  fornula. 

This  probability  adeguately  oodels  those  cases  where 
tenporary  suppression  is  the  only  possible  response  to  a 
given  detonation  of  a round,  it  nay  be  useful  to  expand  the 
model  by  introducing  a second  paired  outcone  at  any  given 
inpacting  round.  Until  now,  we  had  one  paired  outcone: 

Ihe  individual  was  either  suppress,ed 

or  not  suppressed. 

In  addition  to  these  we  consider  a second  paired  outcone: 

The  individual  is  permanently  suppressed  by  being 
wounded  or  killed 

or  not  permanently  suppressed. 

If  we  focus  the  attention  on  modelling  the  effects  under 
this  expansion#  it  is  possible  to  evaluate  a new  probalility 
of  suppression  P(S)  and  a probability  of  kill  P(K)  by  making 
the  following  assumptions; 

Ihe  conditional  probability  of  kill  unlike  the 
condtional  probability  of  suppression  does  not 
depend  on  the  aspect  angle  3. 

It  can  be  represented  by  a smooth  curve  ol!  the 
fcllcving  form: 

l>(K/r)  - exp(-  ~ r^)  (2.20) 

where  a is  a positive  constant  i.e.  H>  0.  This  functicn  was 
selected  on  intuitive  grounds  rather  then  based  on 
real~world  data.  Caking  a vertical  cut  thrcugh  the  surface 
function  above  along  any  angle  one  obtains  the  following 
figure.  Ihe  paraneter  H is  arbitrarily  chosen  as  being 
a»50.  Clearly  this  parameter  is  a function  of  different 
factors  like  terrain  and  aASiUoitlon. 


7/m  3 dlstanoe/meter 


r/n  5 5 10 

for  0,"18O*  foxhole  for  0,«O® 


Figure-4  - FUliCTIOH  P(K/c)  « exp(-  ~ r^) 
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1£  on  the  other  hand  P (K/r)  is  fixed,  0<P(K/r)  <1,  and  the 
equation  is  solved  for  r,  one  will  obtain  concentric  circles 
Hith  iso<-lev6ls  of  probability  of  kill  P (K/r)  around  the 
foxhole* 


The  function  has  the  forn: 
- H*Jln  P(K/r) 


' (2,21) 


The  follouing  fig  5 shows  the.  family  of  functions  for  4 
representative  selective  probabilities  of  suppressions* 

P(K/t)=0*1*i  i=1,2,3,4  (2.22) 

The  parameter  H is  again  chosen  to  be  50* 
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foxhole 


r/a 


Bigura  5 - ISOLIUES  • -H*4a  P(K/r) 


diataace/aater 


a«ii^ef^^riaiaiawi^^ 


Ilte  concentcicity  aioand  the  foxhole  expresses  the 
fact  that  no  natter  in  which  direction  the  person  in  the 
hole  is  looking  the  kill  effect  is  just  deteroined  by  the 
distance « It  would  be  beyond  the  paper  to  verify  this 
assumption  , and  it  wotUd  be  extremely  difficult  to  collect 
data  for  it.  Ihe  reader  must  be  content  for  the  moment  with 
the  earlier  presented  intuitive  argument. 

This  newly  introduced  function  leads  to  a revision 
of  the  probability  of  suppression  P(S)  and  the  evaluation  of 
the  unconditional  probability  of  kill#  P(K).  Ihe  density 
assumed  earlier  for  artillery  hits  is  used  again: 


f(e.r)  . ^ .xp  [-  i (|)2]  ■ 


(2.12) 


Probability  of  suppression  P (S)  but  not  kill  is: 

k 

0»27r  r™® 

P„(S)  = / / P(S/0,r)(l  - P(K/r)]*f(0,r)  r dr  d0  (2.23) 

0«O  r=0 


Bemark:  Ihe  subscript  K at  p^ is  used  to  indicate  that  this 
suppressicn  probability  appears  together  with  the 
probability  cf  kill  P(K). 


Ihe  expression  P (S/9,r)-(1 -P  (K/t) ) means  that 
suppressed  but  did  not  kill  the  individual  at  (r#9) 

Erobability  of  kill  P(K) 


9-2V  r **-®* 

P(K)  - / / P(K/r)*f(0,r)  r dr  d0 

e«0  r-0 


round 


(2.24) 


Ihe  computation  for  yS)  and  ?(K)  runs  along  similar  lines 
as  in  the  earlier  evaluation  of  P(S).  Per  the  specific 
functions  we  find: 
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(2.25) 


0»2ir  r“«»  , - . 1 2 2 

P(K)  » / / exp(->  ^ r ) • y exp(-  -r  (J)  J d(^)  dO 

0-0  r=0  “ 2710^  ^ ° ^ 


Pi/(S)»  / / exp[-  I r^(l  - e cos  0)j  [1  - exp(-  ^ r^)] 

0-0  r-0  ^ ^ 

exp[-  j (^)^3  d(-~)  d0  (2.26) 

2iro 


change  again  the  variahle  of  integration 


2 

r 

2 “ 

z; 

2 

^=0  “ d(^)  = dz 

(2.14) 

P(K) 

6=2 

“ / 
6=0 

/ 2 f-  („  + 2^^1 

2=0  2ito  L “ 0 -1 

(2.27) 

P(K) 

1 

(o^/H)  + 1) 

(2.28) 

P^(S) 


0=2w  , 2 

/ / exp  (1  - c cos  0)  + -^]z  i -i- 


0-0  a-0 


1 


dz  do 


(2.29) 


2iTa 


0=217  z-^«« 


; / 


0=0  z=0  2iro 


2 cos  0)  + I + ds  d0(2.3O) 

^ * 


Becognizing  that  the  first  double  integral  is  exactly  the 
same  as  before  we  can  sioplify  to 


I>^(S) 


'2i 

I 


0«O 


dO 


+ + ^ ) + (-s)  cos  0 

2o  ” ' 


(2.32] 


The  integcatioQ  can  be  done  by  using  the  same  formula  as  for 
the  computation  of  * except  that  this  time  b is:  ' 

(2.33) 


2o‘ 

From  this  it  follows: 


iifter  having  derived  the  desired  probabilities  as  summarized 
on  the  following  figure,  we  are  able  to  evaluate  the 
different  rates  of  temporary  and  permanent  suppression* 

By  assuming  a certain  fire  rate  with  which  the 
artillery  is  firing  in  the  area  where  the  anti-tank  gunners 
ace  located#  the  different  rates  will  have  the  following 
form. 


Tuo-event  model  (unsuppressed-suppressed) 
Bate  of  suppression  1 X.  P(s) 


(2.36] 


?our-event-aiodel  (unsuppressed-supprassed-suivived-Jcillad) 

Bate  of  suppression  X - X.*P,.(S)  (2.37) 

8 t K 


Hate  cf  killing  P(K)  (2.38) 


One  important  rate  must  still  be  developed?  it  is  the  rate 
cf  rise  from  a suppressed  state.  Me  model  this  by  stochastic 
process  methcds. 
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2.  Jate  Of 

The  nodel  so  far  rspresents  a situation  in  which 
there  are  three  possible  outcones.  When  any  independent 
round  iapacts,  the  person  in  the  foxhole  is  assumed  to  he 
either 

suppressed,  - 

net  suppressed,  or 
is  killed. 

Naturally  this  is  only  true  if  we  assume  that  he  is  always 
hack  up  again  when  the  next  round  impacts.  This  fact  brings 
up  the  necessity  of  considering  the  process  along  the  time 
axis. 

Because  of  the  complexity  of  a stochastic  process 
which  reflects  the  behavior  situation  of  an  individual 
exposed  to  arriving  artillery  rounds,  it  seems  useful  to 
start  with  a very  simple  process,  in  which  the  time  of 
suppression  by  a particular  round  is  considered  to  be  fixed. 
In  a next  step,  this  constant  response  time  is  randomized 
over  the  number  of  people  under  consideration,  i.e.  each 
individual  has  his  own  random  time  which  however  is  assumed 
constant  for  the  process  itself. 

Finally,  the  response  time  of  a single  person  say  be 
considered  to  be  a random  variable  coming  from  a certain 
distribution. 

Both  approaches  are  simple,  and  represent  a first  attespt  to 
describe  the  actual  situation  in  different  ways.  Certainly 
the  suppression  time  can  also  be  considered  as  a function  of 
both  above  mentioned  processes  or  of  the  miss  distance  r and 
the  aspect  angle  9,  or  of  a combination  of  all.  But  these 
dependencios  would  rather  complicate  the  mathematical 
derivations  and  laad  beyond  this  study. 

Both  processes  consider  rounds  that  arrive  in  the 


neighborhood  o£  the  foxhole  in  accordance  with  a Poisson 
Process  having  a cate  . The  nodel  consists  o£  two  process 
states.  One  is  the  "up"-state  which  oeans  chat  the  person  in 
the  hole  is  nnsappressed  and  able  to  act  according  to  his 
aission.  Ihe  other  state  is  terned  the  *'dOKn"'’State  which 
results  cron  an  incomaing  round  and  a possible  reaction  to 
it  in  the  fora  of  suppression.  In  this  state  the  perscn  is 
not  able  to  f allfill  his  mission,  he  is  physically  down. 

Ihe  first  process  assuaes  that  the  person  in  the  foxhole 
changes  froii  the  ''‘dOttn'‘-state  to  the  '‘up^-state  only  if  he 
recognizes  a gap  - of  at  least  T time  units  before  the 
appearance  of  the  next  round.  Ihis  means  that  the 
"dcwn"-state  has  a duration  of  exactly  time  T if  no  rounds 
arrive  in  the  time  interval  (0,T).  Ihe  time  T in  connection 
with  such  a process  may  he  called  the  critical  gap. 
[Gaver^p. 48 1 ].  Ihe  following  figure  shows  the  basic 
relationship: 


FIRST  CASS  ; t' > T 


State  0 “upw 
State  1 ’’down** 


T a critical  gap 


f State 


Impact  of  first  Impact  of  second  suppres-sing 
suppressing  >^round  (process  starts  all  ** 
round  over  again) 


T 


0^ 


1 


t'- 


4 


SECOIfD  CASB  : t'  ^ dJ 


A State 


.Impact  of  first 
suppr 
X/''  round 


Impact  of  second  su 


WX  i.d ^ 

^suppressing  round  (process  sta 


over  again) 
T 


■o- 


for  the  tine  being  I is  assuned  fixed.  Houever  it  is  octe 
reasonable  to  describe  it  by  a randon  variable,  since  tbe 
tine  gap  T is  deternined  by  factors  like  learning, 
accustooing  to,  cr  overco;aaing  of,  fear,  stress  etc.  The 
assuaptioQ  of  a fixed  T has  to  be  changed  in  a later  step, 
where  it  will  be  defined  as  a randoa  variable. 

Furtheraore,  the  change  froa  the  ”up«-state  -to  the 
«down"-state  is  accoaplished  by  arriving  rounds  to  which  the 
person  reacts  through  suppression  (‘'down"TState)  with  the 
earlier  computed  probability  of  suppression  P(S)  . 

Now  let  T be  equal  to  the  total  time  of  being 
continually  suppressed.  It  is  possible  tc  define  t in  the 
fclloving  way: 

il  if  no  suppressing  round  arrives  in  (0,T) 
t'  + T'  if  the  next  suppressing  round  falls  (2.39) 


before  T. 


In  this  case  the  random  variable  t'  is  an  independent  version 
cf  the  randoa  variable  x. 

according  to  this  definition  it  is  possible  to  compute  ths 
expected  value  of  x which  represents  the  mean  time  spent  in 
the  ''dowu‘*-state,  or  the  mean  time  of  being  continually 
suppressed.  Ihe  syabol  used  for  it  will  be  ECx]. 

In  order  to  get  the  arrival  rate  cf  suppressing 
rounds  which  contribute  to  the  time  of  being  suppressed,  we 
have  to  multiply  the  fire  rate  with  the  probability  of 
suppression,  F(S) , which  was  evaluated  in  chapter  II.B.I.  He 
obtain 

• ^fVP(S)  (2.36) 
First  we  have  to  state  the  following  two  expectations; 


Elx.lno  round  in  Tj  ■ T 


E[t [round  in  T]  ■ t' 


(2.40) 


i 


ubere  t' » tlie  o£  second  suppressing  round.  Reaoylng  the 
condition  cn  t*  leads. to  expected  tiae  o£  being  suppressed 
under  the  iniluebce  of  incoaaing  rounds. 


EIt] 


EIt] 


-XT  T -X  t*  -XT 

T.c  | e ® X^dt*(t'  -s- E(t1]+BCt1  Cl-e  ®1 


(2.41) 


T»  e + X 


[“X  t 


"j  T 


1) 


-X  T 

+ E(t]  (1  - e ® ) 


Jo 


-XT  . -XT  -XT 

E[t1  - T e ® + f [1  - (1  + X T)  e ® 1 + EIt]  (1  - e ® ) 

A S 


1 


1- 

f 


r 

5 


, .-XT 

E[t]  - T + — ^ (1  - (1  + XgT)e  ® 1 

Xe  ® 

E(t1  - - 1) 


(2.42) 


inserting  « X^*P(S)  we  obtain 


, X P(S)*T 


(2.43) 


I 


Ex-cending  the  idea  for  a fixed  time  gap  by  applying 
it  to  a group  of  people  separately,  ue  are  able  to 
refornulate  the  process  in  the  following  ways 

Assune  a group  of  perscns  in  which  each  individual  sticks  to 
a certain  but  randoi  tine  gap  I uhen  responding  to  a 
suppressive  round,  and  assuae  that  this  randoa  variable  T 
based  on  the  group  is  distributed  with  a density  f(t)  : 

T ~ fp(t)  ■ ^ j (2-44) 
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Ihe  sttDscript  P indicates  the  origin  ot  the  density 
(people).  He  say  express  t which  is  defined  as  before  as 
follows: 


t if  no  suppre^ising  round 
arrives  in  (0,t)  or  t>t 
where  t is  the  tine  of  the 
first  suppressing  round. 
t'+T*  if  t'^t 


(2.45) 


Beoark^T*  has  the  sane  distribution  as  . k sinilar 
derivation  as  before  leads  to  the  following  result. 


~ -i  t*  t nX  t» 

E(T|T=t]  » t / X e dt'  + / (t' +E(t  |T»tJ]X  e ® dt'  ' 

t ® 0 ® 


(2.46) 


”^s“ 

aecall,  the  expression  X„e  is  the  density  of  inconing 

S 

suppressing  rounds. 


"X  t t -X  t*  t 

EItIt-cJ  - tc  t J t'X  e ® dt*  + E(TlT-tJCl  ^ e ® ) 


1 1 "X  ^ 

E(tlT»t]  - ~ ~ e ® + EltlT-tJd  - e ®) 


(2.47) 


solving  for  EC  t|  T=t],  we  find; 
EIt  |T-tJ  ■ 5^  (®  * *“ 


(2.48) 


or  replacing  again  ty  ^3  ■ Xg‘P(a) 


r ••'•q, - 


, ...•TV'Vi 


eIt Ix-t] 


i) 


X.P(S)t 
(e  * 


(2.49) 


So  far  this  value  represents  the  expected  duration  spent  in 
the  "down  state"  given  a particular  tine  gap  t for  a certain 
individual.  Clearly  the  conditional  expectation  has  the 
same  form  as  the  unconditional  expectation  (see  forniula 
(2.43))  for  a fixed  tine  gap  T. 

Beaoving  the  condition  in  equation  (2. 49)  we  find 
the  expected  duration  tine  of  suppression  based  on  the 
considered  population. 


(2.50) 

(2.51) 

00  ^ P(S)t 


E(t]  «•  / E(t  |T»t]»fp(c)dt 
0 


os  1 


For  an  evaluation  of  this  expectation,  the  density  f^(t)  has 
to  be  icnoun. 


h further  variation  of  this  process  leads  to  the 
second  approach.  Here  the  tine  gap  I for  a certain 
individual  varies  randonly  according  to  a particular 
distribution.  This  approach  is  liaited  to  one  individual# 
and  will  not  be  extended  to  a.  group. 
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SiKfpose  an  individual  leacts  vith  a raodoa  tiae  gap  t 

cn  any  incooing  suppressing  round  i«  and  assune  that  this 
randoB  variable  T is  distributed  with 

T-fj(t)  ' (2.53) 

Ihe  subscript  T indicates  the  origin  cf  the  density  (tine) . 
Ihe  duration  t is  defined  as  in  formula  (2.45)  before.  The 
conditional  expectation  can  be  written  as: 

-X  t t -X  t* 

E(t1t  »tl  - te  ® + / (t’  + E[t])  X e ® dt’  (2.54) 

0 ® : 

This  time  t is  conditioned  on  the  individual's  first  chosen 
time  gap  1. 

■ -X  t -X  t 

E[t|T  -t]  - i (1  - e ® ) + E(t](1-  e ® ) (2.55) 

s 

fieooving  the  condition  in  equation  (2.55)  we  can  express  the 
expected  duration  time  of  suppression  based  on  the 
individual's  time  gap  distribution. 

E(t]  ■ / E(t |T,-t]»f_(t)dt  (2.56) 

0 ^ ^ 

E[t]  • f (1  - e ® ) >»•  EttKl  - e ® ) f^(t)dt  (2.57) 


(2.36) 


solving  toe  the  expected  value  SCt]  and  insetting 

Xg  - Aj-P(S) 

we  find: 


e[t1 


-A.P(S)t 

/ (1  - e * ) fj(t)dt 

o 


r -x,p(s)£ 

XjP(S)  [1  - / (1  - e . ) £^(t)dt 


(2.58) 


or  expeessieg  in  Laplace  transform  with  s as  an  argument: 


E[t] 


1 

XjP(S) 


(2.59) 


A further  evaluation  of  this  value  requires  the  distribution 
yt)  of  I. 


comparing  both  expectations, 

Sxpected  duration  time  of  suppression  based  on  the 
population 

(2.52) 


and  expected  duration  time  of  suppression  based  on 
the  individual's  time  gap  distribution 


e(tI 


XfP(S) 


XjP(S)t 


fp(t)dt  - 1 


lie  obsecve  that  they  are  dlffereut.  Houevec  if  we  suppose 

that  both  distributions  f(t)  and  £(t)  are  concentrated  at 

p T 

X«z  (delta  functxon)  it  is  possible  to  reduce  both 
expressions  to  the  very  first  derivation  (2.48)  where  the 
the  tine  gap  I was  fixed. 


Ihe  reciprocals  of  these  expectations  approxisate  the  rate 
at  which  the  foxhole  occupant  (e.  g.  ueober  of  a group  of 
antitank  gunners  returns  from  the  suppressed  state  into  the 
unsuppressed  state)  returns  back  to  continue  his  mission. 


\ 


u 


(2.60) 


where  £Ct]  represents  any  of  the  derived  expectations. 


In  suomary  the  three  rate  coefficients  developed  are; 


Xf‘P(K) 

(2.38) 

Xf-Pj^(S) 

(2.37) 

1 

E[t1 

(2.60) 

If  killing  as  an  additional  event  is  considered,  the  rate  of 
rise  is  computed  with  the  same  formula  (2.60)  except  that 
this  time  ^(S)  is  used  instead  of  P(S) . 


C.  SUPFBESSION  IlOOSI 

The  model  which  will  be  developed  in  this  section  can  be 
set  schematically  in  the  following  framework: 
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I Ibis  analysis  is  restricted  to  the  effects  on  the  anti-tank 
I gunners...  in  this  case  the  X-forces.  it  is  behavior  of 
I . \ nenbers  of  the  X-£orce  that  was  sodeled  previcisly. 

The  approach  to  fornulate  the  situation  of  the  anti-tank 
gunners  is  done  by  using  the  ideas  of  lanchester-type 
eguations  and  their  further  developnent  [ Iayler,1978,p.20]. 

The  differential  eguations  representing  the  nodel  are 
all  deter ninistic  in  the  sense  that  each  of  then  will  always 
yield  the  same  output  for  a given  set  of  input  data.  Even 
though  combat  between  military  forces  is  a complex  random 
process.  These  eguations  shed  light  on  combat  dynamics  and 
may  be  useful  in  defense  planning  studies. 

The  basic  idea  is  that  artillery  forces  use  ''ar6a"-fire 
tc  suppress  or  eliminate  forces  like  anti-tank  gunners, 
"Area**  in  this  context  means  the  fact  that 

the  artillery  unit  "knows"  the  area  in  which  to 
shoot,  but  does  not  know  the  location  of  each 
anti-tank  gunner, 

the  anti-tank  gunners  are  "invisible"  to  the 
artillery  unit. 

If  we  further  assume  homogeneous  forces  of  X,  it  is  possible 
to  set  up  differential  eguations  which  model  the  rate  of 
change  of  the  X-forces: 

dX  (t) 

-ft 

dX  (t) 

-dT-  * 


aeiark:  Claacly  A bare  is  takes  £koa  the  U-event-aodei.  The 


variables  and  X^Ct)  are  respectively  the  nuiaber  of  the 

X'-forces  which  ace  either  active  in  the  foxhole  (able  to  use 

guns  and  anti-'tank-veapons) , }^(t) , or  suppressed  , X^Ct). 

In  this  simplified  structure#  the  two  equations  describe 

the  most  essential  factors  of  the  assumed  situation.  They 

are  mathematically  approximate  , because  the  solution  of  the 

two  differettial  equations  will  furnish  the  num'ber  of 

gunners  active  on  the  battlefield  and  the  number  of  gunners 

supppressed  at  any  point  in  time.  The  derivation  uses  the 

¥ 

Laplace  transformation: 


^ {X(t)}  = x(s) 


(2.63) 


we  perform  Laplace  tcansformaticn  upon  (2.61)  and  (2.62) 


sx^(s)  - X^(0)  " (Xg  + Xj^)  x^(s) 

+ X X 
u s 

(s) 

(2.64) 

sx^(s)  - X^(0)  = \x^(.s)  - A^x^(s) 

(2.65) 

Solving  these  two  equations 

for 

X(t) 

a 

and  X (t)  and 
s 

translating  them  back  to  XJt)  and 

a 

time  dependent  quantities: 

Xjt) 

gives 

the  desired 

-X  (0)  » - (X  + X,  + s)  X (s)  + X X (s) 
a s K a u s 

(2.66) 

-X  (0)  » X^x  (3)  - (X  + s)  X (s) 
s 8 a u s 

(2.67) 

-X  (0)  [X  + s]  ■ -(X  3)  + X (X  + s)  X (s) 

duus 

-X^CO)  X^  - XsV^(s)  - (X^  s)X^  Xg(s)- 

-X^(0)(X^  -f  s?  - X^(0)V^  - XgX^x^(s)  . (X^  + Xj^  + s)(X^  + s)  x^(a) 

X (0)  X 4-  X (0)  X + sX  (0) 

(X  + X.  + sHTTsT^nTx 

S Ic  U 8 U 


(2.68) 


BUltiplyiltg 
(X^  + X.  + s) 


(2.  66}  and  (2.67) 
respectively  wo  receive: 


«3(s) 


X (0)X  + X (0)(X  + Xj  + sX  (0) 
a s s s K 8 

(X  +-X.  + s)(X  + s)  - XX 
« k u s u 


with  Xg  and 


(2.69) 


Suppose  that  at  t=0  the  number  of  active  gunners  (able  to 
uatch  and  shoot)  is  equal  to  X^and  the  number  of  gunners 
being  suppressed  (down  in  the  foxhole)  is  equal  to  0^ 
for  i.e.  t-0 


X^(0)  = Xq 

X (0)  » 0 
s 

the  equations  for  X(e}  and  X (s)  can  be  rewritten: 

a s 


(2.70) 


ii^(a) 


0 u 


T 

+ s(X  + X,  + X ) + X,  X s + s(X  + X.  + X ) + X.  X 
sku  ku  sku  ku 


X.X 

X (s)  = ^ (2.72) 

® 3+s(X+X,  + X)  + X,  X 

s k u k u 

Using  the  Laplace  correspondence: 


1 

— ^ ^ 

s + (-a-b)s  + ab 


for  a b 


at  , be  ) 
r ufi  ~ be  I 

oC  a - b J 


2 

s 


3 

+ (-a-b)s  + ab 


(2.73) 


for  a b 


Hhece 


a • b a A,  X 
® k u 


and  -a  - b = X + x,  + X 
s k u 


(2.74) 


bence  a » f (-(x^  + + x^) ) + /(X^  + x^^  + x^)‘  - AXj^X^ 


and  b ” ■*■  ^ ^ 


a k u 


k'u 


(2.75) 


Ibe  equation  for  X (t)  and  X (t)  than  are: 

a s 


- (X  + b)e^‘^] 

a a "*  0 U U 


(2.76) 


X 

w /j.  \ _ 0 . f _ at  be  1 

X^(t)  = r X^[e  - e J 

s a “ D s 


(2.77) 


Since  X^,  and  X^  have  tc  be  always  greater  or  equal  to 


zero 


s k 

X > 0 
s — 


(2.78) 


Ihe  constants  a and  b are  always  negative  and  real  numbers, 
a i 0 

(2.79) 

b £ 0 

Ibis  leads  to  the  basic  shape  of  the  functico 


X^(t)  - f(c) 
d 


Xg(t)  - f(t) 


(2.80) 


Shown  in  the  following  figure. 
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In  discussing  the  functions  X (t)  and  X(t)  the  folloving 

A 8 

pcoperties  can  be  seen* 


for  t » 0 


x^(0)  « -^-3^  [(X,^  + a)G®*°  - (Xy  + b)e^*°l  - (2.81) 


,,  fn\  0 , r a*0  b*0, 

— K •*  ® 3=0 

8 a - D s 


(2.82) 


for  t » » 


X («)  = lira  ^ [(X  + a)*e^*' - (X  + b)*e^*']  = 0 (2.83) 

a a - b u • u 


^ . 0 . r at  bt, 

X (“»)  = lira  - — -~r  X [e  - e j = 0 
s ^^.^a-b  s 


(2.84) 


for  t St  (i.e.oax  values  for  X (t)  and  X (t) ) 
sax  a s 


X (t)  has  no  max  in  (0,»),  i.e.  max  value  is  at  t = 0 
a 


*a'‘>  • 


(2.85) 


X (t)  has  a max  at 
s 


in(a/b) 

t ® with 

max  b - a 


(2.86) 


Xrt  at  ^ bt 

X (t  ) » ^ X [e  - c < X (2.87) 

s max  a - b s 0 


Ibe  foregoing  aodeX  assumed  a rate  of  killing  in  addition  to 
a suppression  rate  and  a rising  rate.  It  is  possible  to 
simplify  this  model  by  leaving  out  the  third  rate.  It  will 
be  shown  in  chapter  III.  that  this  rate 
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for  such  a scenario  aay  be  very  soaXl  in  ccaparlson  to  the 
other  tvo  cates. 

For  this  reason  we  can  state  that  leaving  out  the 
killing  rate  vill  not  drastically  oversiaplify  the  BOdel, 
yet  it  uill  sinplify  the  ccnputaticnal  procedure  in 
obtaining  the  wanted  dependency  between  tioe  and  the  nunber 
of  gunners  active  or  suppressed  on  the  field.  ’ 

This  leads  to  the  following  set-up  of  differential 
eguations: 


dx  (t) 
a 

dt 


s a 


+ X X (c) 
us 


(2.88) 


dX  (t) 
s 

dt 


X X (t)  - X X (t) 
s a us 


(2.89) 


where  the  total  sum  of  people  either  suppressed  or  active  on 

the  battle  field  is  egual  to  a constant  1 (no  killing) 

0 

i.a.  X (t)+X  (t) 
o a s 


which  enables  us  to  rewrite  the  equations! 


dX  (t) 

-f-—  » -X  X (c)  + X (X_  - X (t)) 
dt  s a u 0 a 


(2.90) 


dX  (t) 

-fi-  - 


(2.91.) 


Using  again  iaglacn  transfornation  as  earliex; 

X X 

sX^(s)  - Xq(0)  ■ (2.92) 

aX^(a)  - X^(0)  - tA  . - i„X,(a)  (2.93) 


I 


S6 


X (s) 
a 


X (s) 
s 


sls+Uj  + X|j)l  s + i\^+XJ 


>A  , 

sis  + U.  + 1, .))  s + (».  + »„) 

S u 9 u 


using  the  Lagiace  transfors: 


=7T7 


(2.94) 


for  both  equations  we  find  X^(t)  and  X^t)  fct  the  simplified 
model. 

The  rate  has  to  be  derived  this  time  from  the 

9 

pcobability  of  suppression  P (S)  without  Killing  included  as 
a possible  event.  We  recall j 


, (2.19)  X » X,  • P(S)  (2.36: 

S £ 


o 2 B , ^ K 2 

2ira  I 1 + — s-  - c 

' 20'^ 


hence 


X (t)  » X^Xq  ^ (1  - exp[-(Xg  + Xy)t])  + X^(0)  exp[-.(Xg+X^)t] 

^ ^ (2.95: 


• ^^0  IT+T’'  ‘ exp[-(Xg  + X^)u] 

^ . n ® (2.96: 

for  t = 0 

X (0)  » X,.,  X^(0)  = 0 (2.97 

a s 


The  eguaticn  for  X (t)  and  X (t)  then  are: 

a s 


s u 


(2.98 


X.to  • Xs  T 


J.  \ 


(1  ~ axpHX  + 


(2.99 


Beaark:  foe  t the  sus  of  both  assyn^totic  values  of 

X(t)  and  X(t)  is  equal  to  x^. 
a s 0 

Ihe  basic  shape  of  these  tvo  functions  X(t)  and  X(t)  can  be 

d 3 

seen  in  the  following  Fig  10. 

If  ve  cenparo  this  result  to  the  result  on  figure  9 we 
can  see  that  the  addition  of  killing  to  the  model  has  an 
influence  On  the  shape  of  the  functions.  So  is  e.  g.  the 
value  for  X(t)  and  X(t)  in  the  first  model  (figure  9)  for 

3 S 

large  times  approximately  zero,  while  in  this  case  here 

(figure  10)  , the  total  number  of  people  on  the  field  (XJt)  ♦ 

X(t))  is  always  constant,  i.  e.  both  functiens  approaches  to 
s 

a limit  value  not  equal  to  zero. 


III.  FITIINS  EXPERIMEHTAL  DATA  DESCRI|INg  SOPPRE^SIOH 


A.  £XP£BIM£NIAL  DESIGN 


general  Aspects 

Ihe  natheoatical  oodel  developed  in  the  chapter  before 
shall  now  be  supported  by  an  experiment  vhich  was  conducted 
by  the  US  Army  Combat  Developments  Experimentation  Command 
(CDEC)  # Ecrt  Ord,  California. 

The  experiment  described  here  belongs  to  a series  of 
similar  experiments^  all  dealing  with  the  objective  of 
ccllecting  and  analyzing  data  of  the  suppression  process. 
The  data  of  this  particular  part  of  the  CDEC  experiment  were 
based  on  the  relationship  between  suppression  and  distance. 
This  analysis  tries  to  make  use  of  the  data  by  analyzing  the 
relationship  among  suppressionr  distance,  and  angle.  In  a 
further  experiment  conducted  by  CDEC,  the  objective  was  also 
to  include  the  angle  as  an  additional  variable  for  the 
suppressicn  effect. 

2.  setup  Realisation 

The  data  were  taken  from  a part  of  the  experiment  which  was 
executed  at  Fort  Hunter  Liggett,  California. 

four  foxhole-bunkers  were  constructed  which 
guaranteed  the  safety  of  the  players  as  well  as  reproduced 
the  real  scenario  as  close  as  possible.  Their  tops  were 
below  ground  level,  and  covered  with  several  layers  of  wire 
■esh  and  steel  plates.  This  provided  overhead  protection 
from  fragmentation.  Each  bunker  was  equipped  with  a mirror 
system  which  allowed  the  player  to  view  events  in  front  of 
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bis  position. 

A pop-up  silhouette  vas  installed  fcrwaid  o£  each 
player  position.  The  player  vas  able  tc  control  the 
silhouette  as  veil  as  the  cover  of  the  eirroc  systea,  i.e. 
uhen  the  cover  vas  opened  (allovi::^  the  player  to  viev 
dovn-range) , the  silhouette  vas  in  exposed  posture.  He  vas 
asked  to  respond  as  if  he  vould  be  in  the  position  of  the 
pop^'Up  silhouette. 

Each  bunker  vas  connected  to  controll  bunkers  by 
comaunication  and  instrumentation  vires  and  power  bunkers 
for  data  recording  and  supply.  In  the  forefield  of  the 
bunkers^  different  types  of  projectiles  or  equivalent 
charges  (31  mm, 105  mm, and  155  na  aaong  others)  were  placed 
aud  randomly  detonated  one  at  a time  with  the  time  between 
detcsations  randomly  selected  from  three  possibilities  of 
ten,  fifteen,  and  twenty  seconds.  The  figure  on  the  next 
page  shows  the  schematic  setup  of  the  rounds,  the  location 
of  the  bunkers  and  the  angle  of  sight  for  each  foxhole.  it 
can  be  seen  that  the  explosions  were  visible  to  some  players 
hut  not  to  others.  Since  each  type  of  ammunition  .has  a 
different  lethal  radius,  it  was  necessary  to  have  different 
range  configurations  for  each  type.  The  aspect  angle  and 
the  miss  distances  from  the  foxholes  to  the  different  points 
of  detonation  are  summarized  in  appendix  A for  each 
ammunition  separately.  Since  all  of  the  projectiles  used  in 
this  part  of  the  experiment  were  statically  detonated,  it 
was  not  possible  to  model  the  kinetic  contribution  to  the 
terminal  effects.  In  order  to  keep  the  fragmentation  pattern 
as  close  as  possible  to  i;hose  of  incoming  rounds,  the 
projectiles  were  supported  with  sandbags  [ 3uppex,p«A'*12]. 


m 


Figure  11  - SCHFmTIC  RAJIGE  FOR  01/105/155  mm 
SUPPRESSION  EJCPaUMENT  (SUPPBX  II) 
CDEC 
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lateral 

dlstaace/m 


The  players  were  divided  into  two  foufoan  teaas* 
For  each  trial#^  aeobers  of  one  teao  occupied  individual 
foxholes  and  provided  the  perforaance  data  for  that  trial. 
The  aission  of  the  players  was  to  track  moving  target  tanks 
by  operating  the  periscopes.  This  aission  was  interrupted  by 
the  players  responses  to  detonations  in  case  he  was 
suppressed  (change  of  his  state  to  "down") . The  state 
change  and  the  pericd  of  suppression  were  automatically 
recorded. 

3.  presegtaticn  0|  The  Da^a 

Appendix  A summarizes  the  data  which  are  the  basis  for  the 
succeeding  data  analysis  and  serve  tc  evaluate  the 
paraaeters  and  rates  for  the  model. 


E.  DATA  ANALISIS 

“J*  Parameters  21  Ihe  Model 

The  data  presented  in  the  section  above  are  separated  into 
three  different  sets: 

Set  one  consists  of  data  related  to  the  conditional 
frequency  of  suppression  #(S/9,r)  and  the  time  of 
suppression  when  rounds  of  caliber  81  mm  were  fired. 
Specifically,  if  n trials  were  made  under  conditions 
(r,8) , the  number  s^cf  suppression  was  tabulated. 
Then  ^/n  is  an  estimate  of  P(S/r^8j,). 

Sets  two  and  three  provide  the  sane  data  except  that 
they  are  related  to  105  amm  rounds  and  155  an 
rounds. 

The  aain  part  of  this  analysis  is  to  make  use  of  the  data  in 
order  to  estiaate  the  paraaeters  of  the  aodel  described  in 
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chapter  II,  and  to  see  hot?  well  the  real  world  situation  can 
be  described  by  the  nodel  developed  earlier. 

The  model  is  fitted  (numerical  values  of  the 
parameters  are  determined)  by  the  method  of  least  squares. 
[Oe  Groot,  p.S27'*538j.  Nith  the  basic  model  being  of  the 
form: 


®i“i  * Sa'i 


(3.1) 


And  the  data  are  analyzed  on  the  following  two  different 
assumptions;  the  second  is  certainly  the  most  realistic: 

Homoscedasticity  of  the  data,  i.  e.  the  error 
variance  of  Err  is  assumed  to  be  constant. 

Hetexoscedasticity  of  the  data,  i.  e.  we  will  assume 
that  each  error  term  Err  is  distributed  with 
variance  , where  the  variance  is  not  constant  over 
observations.  Errors  are  also  assumed  tc  be 
independent.  If  convenient,  error  terms  will  be 
assured  to  be  approximately  normally  distributed. 

These  assumptions  imply  different  regression 
methods.  For  the  first  assumpticn,  the  result  of  unvf^ighted 
single  step  and  iterative  regression  methods  will  be 
presented.  On  the  assumpticn  of  heteroscedasticity , two 
special  iterative  methods  will  be  used.  When  we  allcw  for 
heteroscedasticity,  ordinary  least-squares  estimation  places 
the  same  weight  on  the  observations  which  have  small  error 
variances  as  on  those  with  large  error  variances.  By 
applying  a weighting  regression,  it  is  possible  to  adjust 
for  the  heteroscedasticity.  so  the  two  announced  methods  for 
the  second  assumption  are  iteratively  weighted  least  squares 
regression  methods. 
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In  order  to  prepare  the  fornula  of  the  aodel  for  the 
regression#  we  apply  a log*‘transfornation  to  the  equation 
(2.1)  in  chapter  II  B.1.  This  is  a convenience#  for  it 
transforos  the  problem  to  one  of  linear  fitting,  we  obtain 
the  following  result: 


P(S/9,r)  = exp[-  ^ r^(l  - e cos  0)  3 (2.1) 

in  P(S/9,r)  * - ^ ^ cos  0 (3.2) 

Using  the  notation 


y = In  P(3/0#r) 

u a r^  (3.3) 

V =r^-cos0 

the  equation  can  be  applied  to  each  datapoint  1 and  the 
equation  can  be  rewritten  as; 


K K 


(3.4) 


The  model  can  then  oe  expressed  as  follows: 


j 

> 


where  the  uukncwn  parameters  are:  6^ 


(3.1) 
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and  the  randcn  variable 


I 


fV 


> 

1 


1'' 


't 

i; 

f 


1 


1 


- in  P(S/0^,t^)  (3.6) 

Ihe  “hat''  on  the  letter  P represents  an  estiaate  of  the 
probability  of  suppression  at  (r,9).  the  fitted  regression 
plane  for  this  function  has  the  fora:  ''' 

y » + ^2^  (3^7) 


The  log-transformation  of  the  conditional  frequency 
of  suppression  ^ (S/e,r}  causes  difficulties  because  of 
experimental  results,  which  lead  to  an  observed  frequency  of 
zero  suppression.  This  fact  influences  the  estimate 
i(S/0,r),  that  is  used. 

Having  specified  the  model  it  is  possible  to  apply 
the  following  different  regression  methods  on  the  assumption 
of  homoscedasticity: 

In  other  words,  suppose  n^^  observations  were  made  at 
experimental  conditions  and  of  these  s^  were  successes. 
Then  first  consider  the  estimate 

- in  ks/r^,0^)  - in(s^/n^)  (3-8) 

where  s^  is  the  number  of  suppressions,  and  Oj^  is  the  number 
of  exposures  to  suppression,  under  condition  1. 

Consequently,  the  transformed  response  to  conditions  (2^,^) 
may  be  In  (B/n^^)  a - «>  formally,  causing  eobarassment.  ye 
will  subsequently  suggest  alternative  ways  of  dealing  with 
this  problem. 
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flatbed  1:  (z€i:o>probabilities  ace  oaitted) 


With  this  aetbod  , the  datapoints  which  had  a 
fraguency  of  suppression  P (S/r,6)  =s^/n^  -0  are 
deleted.  Consequently  the  following  data  points  of 
the  original  sample  of  appendix  A as  shown  in'  figure 
fig  12  were  net  considered. 
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Uethod  2 and  3 (Clustering  of  data  points) 


One  or  nore  zero  data  points  i , i.e.  P(S/ei>r)  s o 
i=1,2/...a  will  be. clustered  with  one  other  data 
pcint  j which  has  a nonzero  conditional  frequency  of 
suppression^  i.e;  P(S/9,r)  # 0 and  which  possesses 
the  closest  distance  and  angle  6^  to  the  zerc  data 
pcint  (s)  • 


By  clustering  (o+l)  data  points  together  and  taking 
the  mean  of  the  frequencies  , e and  r within  each 
cluster,  we  are  able  to  replace  the  data  points  of 
the  cluster  sets.  This  procedure  is  inevitably 
somewhat  arbitrary. 


^ . nH-X  ^ 

P(S/0,r)  “ ^ I P^(S/0,r) 
i»l 

1 1 

» ■ ST  , Vi  ' ‘ ^ ill  'i 

1=1 


(3.9) 


The  bars  above  the  symbols  represent  the  cluster 
average.  with  these  averages  two  different 
replacement  procedures  can  be  applied: 

Bach  cluster  will  be  replaced  by  its  cluster 
average.  This  reduces  the  data  sample. 

Bach  cluster  point  will  be  replaced  by  its  cluster 
average.  Consequently  the  number  of  data  paints 
stays  unchanged. 

figure  13  lists  the  replaced  data  points  of  the  samples  of 
appendix  h. 


69 


•n  tn  «n 
?*i  {M  (M  rj 

r«.  r- 
r>  w esj 
rH  rH  H 

d d d d 

fH  rH  rH  tH 
H «H  rS 

ir»  m «r\ 
O'  CN  0\ 

lOA 

104 

tn 

r-i  tH 

OOO 
cn  CO  CO 

r-l  r-<  H 

M C'4  CM  N 
OOOO 
rH  rH  iH  H 

200 

200 

Co 

OOO 
CM  CS  CS 

125 

125 

ohOvOvO  OOOO  OOO  yo  \D  •j-sf'SfsC  OO  CMCSC'J  mco 

Irtcnn  V0\0«0  r-4r-l  MCO  SfsT'S'  V0^0^<»^0  i-IHr-t  MN 


127 

110.25 

95 

104 

'O 

CO 

iH 

130 

102 

200 

207 

125 

tr\ 

yO 

o 

d 

'O 

>0 

o 

CM 

fO 

cn 

sj- 

yO 

CO 

«o 

VO 

iH 

<n 

CO 

CO 

yo 

VO 

sC 

cn 

O' 

CO 

fs. 

O 

o 

O 

rH 

o 

o 

o 

iH 

o 

• 

• 

• 

♦ 

• 

• 

* 

» 

• 

• 

o 

o 

O 

O 

O 

o 

O 

o 

o 

o 

IN. 

00  CO  IN.  00 

in  lo  m 

00  O 

iri  in  o 

eo  00  IN.  r>. 

o o 

CM  C^  O' 

m m 

<N  M CS 

O O <-(  o 

O'  O'  c^ 

o o 

IT)  l-l 

n n 

O O O'  c^ 

o o 

OOO 

CM  CM 

«-<  H cH 

iH  H (H  iH 

f-l  rt 

H iH 

^ iH  fH 

rt  rM 

• M CM 

CM  CM  CM 

fH  H 

uo  m uo 

• « t 

\0  \C  yO 

vO  vO  rH  IN. 

OOO 

fV 

vo  r*» 

VO  MJ  O' 

vD  vO  rM  H 

o O 

fv. 

CO  CO 

CO  CO  CO 

nS-  -C  <■  CM 

VO  \0  VO 

^ fH 

CO  CO 

“a-  >»  sr 

m m fN, 

H 

CM  CM 

CO 

CO 

<n 

c*. 

CO 

e\ 

00 

<r 

iH 

H 

ir» 

H 

H 

H 

fH 

CO 

H 

* 

• 

• 

« 

« 

* 

» 

• 

• 

• 

OOO 

OOOO 

OOO 

O O 

o o 

OOO 

o <r>  o o 

o o 

OOO 

o o 

CO  o \o 

m CM  O'  'fl’ 

Mf  *H  OJ 

•O'  »o 

O'  o 

CO  sf  m 

CO  O'  O >H 

CO 

c^  m 00 

MT  Cn 

CO  CO  CO 

CO  CO  CM  CM 

CO  CO  CM 

iH  H 

H CM 

CM  CM  CM 

CM  CM  CO  CO 

^ iH 

H H 

i 

tH 

1 

try 

I 

«n 

CO 

O 

in 

K 

H 

Pigare  13  - caANGEO  oaia  points 


70 


flethqd  4 (data  transforoatioa  by  Cox) 


This  method  Is  analogous  to  one  suggested  hy  D.  B. 
Ccx  [Cox#  p.333.  The  derivation  for  the  variable  y^ 
shewn  on  these  pages  can  be  applied  correspojidingly 
in  the  following  manner: 


(3.10) 


where  the  random  number  is  the  number  of 
suppressions  (successes#  given  n^^  trials) . 

The  constant  a#  which  represents  an  unbdasing 
adjustment#  is  derived  in  appendix  B.  It  is  taken  to 
be: 


a • Y (1  - ks/0,r))  . 


(3.11) 


This  is  suggested  by  an  auxiliary  analysis  similar 
to  that  of  Cox* 

Hence  the  transformation  results  in  tha  formula: 


h " 


+ I (1  - P^(S/83r))^ 


(3,12) 


This  formula  enables  one  to  include  data  points  that 
involve  zero  observed  freguencies  of  suppression; 
the  eabarassaent  of  taking  the  logarithm  of  zero  is 
no  longer  present. 

n 


- a 


Method  3 (unweighted  itecative  regressiqn) 


As  an  extention  of  the  nethod  discussed  last#  this 
aethod  replaces  the  value  of  P(S/e#c)  in  the  focoula 


. At  4 (1  - ?t(S/e.r))j 


(3.13) 


After  each  iteration  by  an  estimate  of  the 
probability  of  suppression  P (S/e#r)  # the  number  of 
iterations  vas  determined  by  the  appearance  of 
convergence  of  the  parameters  X and  e. 


The  two  following  methods  are  weighted  regressions  which  are 
necessitated  by  the  earlier  described  neteroscedasticity. 


Method  6 (weighted  iterative  regression  Var  {Err)~r^) 


Ihis  method  is  suggested  by  the  fact  that  the 
observed  variability  in  residuals  increases  with  r. 
[Pindyck#  p.  100]  Since  the  error  variance  is  not 
known#  it  is  reasonable  to  assume  the  existence  of  a 
simple  relationship  between  the  error  variances 
Var  (Srr^,)  and  the  values  of  one  of  the  explanatory 
variables  in  the  regression  model*  In  this 
analysis#  the  distance  from  the  foxhole  to  the 
explosion  was  chosen  as  the  important  explanatory 
variable.  By  using 


VarlErr^l  ~ 


(3.14) 


the  resulting  regression  foroula  is: 


(VarlErrJ) 


1/2 


^ (Var(Err^])^^^ 


+ e. 


(Var(Err^l)^^^ 


(3.15) 


Err^ 


(VartErr^l) 


172 


This  equation  can  he  reduced  to  the  folloving 
regression  function: 


y'. 


^1^ 


(3.16) 


where 


in 


Sj  + i(l  - 


*1 


) 


(3.17) 


Afterwards#  the  iteration  procedure  is  equivalent  to 
netbod  5.  This  method  tends  to  sake  the  variance  of 
the  residuals  around  the  fitted  line  of  aore  nearly 
constant  variance;  estimates  of  the  aodel  paraaeters 
should  be  thereby  iaproved. 


tiathod  7 (weighted  iterative  regression) 
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Using  the  error  variance  derived  in 

which  has  the  fora: 

appendix 

1 - P (S/0,r) 

Var[in  f^l  “ 

^ P^(S/e.r) 

(3,18 

fhere 

> 

(3.19 

the  regression  formula  will  be  modified 

as  follows: 

(Var(Jln  f^J) 


where  y is  again: 


(Var(in  f^]) 


/S.  + i (1  - P,(S/0,r))\ 

^ } 


(Var[ln 


(3.20 


(3.12 


as  derived  in  nethod  4.  ihe  resulting  equation  has 
the  form;  “ 


S,  + 4 (1  >•  P,(S/9,r)) 


1 - P,(S/0,r) 


n^*  P^(S/0,r) 


^/l  - P,(S/9,r)V^^ 


n,  • P,(S/9,r) 


‘ /I  •*  P^(S/9,r)\ 
\n^»  P^(S/0,r)/ 


The  iterative  procedure  in  this  oetbcd  is  perforoed 
hy  using  the  estiiate  of  the  conditional  probability 
of  suppression  as  input  for  each  succeeding 
iteration.  The  nuaber  of  iterations  was  deteroined 
as  in  aethod  5 and  6 by  the  apparent  convergence  of 


the  estinatocs  of  the  paraoetecs  K and  e • 


In  order  to  perfora  the  multiple  linear  regression,  ve 
switch  over  to  matrix  notaticn,  where  we  can  write  the 
normal  equations  in  the  follcuing  form:  ^ . 


fp  A . m 

z“*z*e  » z ‘Y 


(3.22) 


Ihe  matrix  Z is  the  design  matrix  which  consists  of  vector 
0,  the  square  of  the  distances  and  V,  rhe  product  of  the 
cos6  and  the  square  of  the  distances.  The  variables  y and  B 
represent  a {sample  size  l^l)  and  a (1x2)  vector  respectively, 
for  the  experiments  analyzed,  the  entries  of  the  design 

A 

matrix  were  obtained  by  solving  the  normal  equations  for  B* 
He  find: 


• T — 1 T 
B « (Z  *2)  ^‘Z  ‘Y 


(3.23) 


Bemark:  Capital  letters  used  for  matrix  notation. 

Ihe  parameters  K and  e can  be  evaluated  by  using  the 
equations: 


A 


K - 


£ - Bj-K 


(3.24) 


AS  a modification  of  the  supprassicc  function  (2.1)  * ' 
in  chapter  II.  B.  1 the  function  (2.4)  in  the  same  chapter 
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was  cons i dec ed 


h • 


P(S/0,r)  « L*exp  ~ r^(l  - e cos 


(2.4) 


The  iaclusion  of  the  parameter  L allows  for  the  possibility 
that  suppression  may  not  occac  for  shots  that  fall  in  very 
close  proximity  to  the  foxhole. 

dethod  d (sisilar  to  method  4 but  using  model  above)- 


In  this  method  the  regression  procedure  and  the 
method  4 were  applied  to: 


(3.25) 


where 


■=•  in 


+ I (1  - P^(S/9,r))^ 

\ / 


cos  6 


®l--t 


(3.3:6) 


®2“t 


83-0* 


The  scrmal  eguatlons  have  tne  same  fora  as  befcce: 


B - (2^*Z)''^«2^*Y 


(3.23) 
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but  It  the  design  matrix  consists  this  tine  of  the 
vectors  (J«V,  and  a colunn  vector  consisting  only  of 
ones.  The  resulting  0 vector  consists  now  of  3 
elenents,  fron  which  it  is  possible  to  evaluate  the 
different  paraseters  K,  s » and  L. 

The  calculation  of  the  parameters  for  all  methods  was 
performed  fcr  each  type  of  ammunition  separately  by  applying 
the  computer  language  API.  The  programs  which  where  used 
for  this  regression  are  listed  in  appendix  X.  The  appendixes 
b through  G show  the  analysis  of  the  residuals  after  the 
regression  ci  the  original  data  of  appendix  A. 

Ihe  analysis  for  each  ammunition  type  was  performed 
according  to  the  same  pattern,  and  includes  the  following 
steps. 

Por  each  earlier  developed  method  (1  through  7)  : 
Determination  of  beta(1}  and  B£XA(2) 

Determination  of  K 
Determination  of  e 

Plot  of  the  residuals  as  a function  of  y as  defined 
in  each  method  earlier.  The  APL  function  which  was 
used  has  the  name  SCAT  and  belongs  to  the  library 
package  OA  2 3660  (available  at  B.  Church 
Computer  Center). 

Plot  of  the  residuals  of  the  regression  as  a single 
array  with  the  function  BOXfLOT  of  the  same  library 
package^  The  plot  characterizes  the  guartiles,  the 
interquartile  distance,  the  median,  data  points 
inside  and  outside  the  1 and  1.5  interquartile 
distance  and  outliers  [hcNeil,  p.13  and  71,72]. 

huaerical  values  of  the  residuals. 
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Histograns  £oc  the  residuals  of  the  regression  which 
show  the  relative  freguency  and  statistical 
features. 

For  method  8,  which  is  an  analysis  based  on  different 
assumptions  for  the  probability  cf  suppression  P(S/8,r},  the 
same  analysis  was  performed,  with  the  exception'of  the 
histograms  fcr  residuals. 

Ihe  methods  1-7  describe  the  attempt,  to  master  the 
problem  of  zero-probabilities  in  the  original  data  set  as 
well  as  the  problem  of  a possible  heteroscedastisity. 

For  the  later  methods,  in  which  more  appropriate 
statistical  tools  were  used,  the  systematic  structure  of  tne 
residuals  seems  to  disappear  up  to  a certain  point.  The 
residuals  concentrate  themselves  more  and  more  symmetrically 
around  their  mean  and  median  (compare  particular  the 
boxplots  and  the  histograms  in  the  appendixes  D,  E,  F,  and 
G)  . 

The  relative  large  remaining  range  of  the  residuals 
is  determined  by  single  outliers.  The  analysis  showed 
differences  for  different  kinds  of  ammunition.  Among  the 
three  ammunitions  considered,  the  analysis  of  81  mm  showed 
the  smallest  spread  for  the  residuals,  for  almost  all 
regression  methods. 

The  appendix  H shows  for  the  iterative  regression 
methods  5,6,  and  7 the  development  of  K and  e for  the 
different  kinds  of  ammunition,  in  each  method  , a 
convergence  of  the  values  K and  e with  increasing 
iteration  can  be  observed.  The  starting  value  for  K is  in 
method  5 smaller  and  in  method  6 and  7 larger  than  its 
corresponding  value  after  convergence  is  obtained.  This  is 
true  for  all  types  of  amoiunition.  An  equivalent  observation 
can  be  made  for  the  e -values. 

A graph  for  the  iterating  e for  method  5,  155  an 
could  not  be  plotted  because  of  the  very  small  change  of  the 
£ -values  along  each  iteration. 
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The  third  page  of  appendix  H displays  the  numerical 
values  for  K and  e produced  in  each  method  for  the  three 
different  types  of  ammunition. 

The  analysis  of  the  parameters  K and  ^ showed  the 
following  result; 

for  the  iterative  methods  (!j,6»7)  the  K and 
Srvalues  approached  with  increasing  iteration  a limit  value 
(see  page  one  and  two  of  appendix  H)  ; this  occured  after 
about  the  sixth  iteration.  The  final  values  for  these 
aethcds  are  also  presented  in  the  tableau  of  appendix  fl 
(third  page) . The  plots  for  K and  e show  different  shapes 
for  different  methods  and  partly  also  for  different  types  of 
ammunitions  , The  approach  to  the  final  value  may  occur  from 
a relative  snail  or  a relative  high  value.  Among  the  three 
ammunition  types  considered,  the  analysis  of  81  mm  showed 
tue  smallest  spread  for  the  residuals,  for  almost  all 
regression  methods.  The  scale  parameter  K influences  the 
probability  of  suppression  P(S/9,r)  as  stated  under  formula 
(2.1)  in  chapter  in  the  following  manner; 

An  increasing  value  c*f  K leads  to  an  increase  of  the 
probablility  of  suppression  P(S/0,r)  .It  is  reasonable  to 
assume  that  the  suppression  effect  increases  with  the 
increase  of  the  caliber.  This  behavior  was  confirmed  by  the 
data  analysis;  except  for  method  one  and  three  where  an 
inversion  could  be  observed  between  the  K* values  for  81  mm, 
105  mm,  and  155  mm  and  for  81  mm  and  105  mm  respectively. 

This  distortion  results  from  the  fact  that  in  irethod 
one  all  data  points  with  probabilities  of  suppression  equal 
to  zero  are  disregarded  and  in  method  three  the  cluster 
procedure  emphasizes  the  average  values  produced  by  the 
clustering.  He  prefer  the  latter  methods. 

Contrary  to  these  observations  on  the  K*“  values,  a 
general  trend  for  the  e>values  can  not  be  related  to 
different  methods  or  different  hinds  of  ammunition.  The 
most  reliable  value  of  K and  e seems  to  be  found  by  method 
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7.  This  oetbod  considers  the  different  variabilities  for 
different  sets  of  distance  and  angle. 

In  the  course  of  this  paper,  the  value  of  K and  £ 
derived  in  sethod  7 are  taken  as  input  for  developing  the 
rates  for  the  model  of  section  Il.B. 

Ifaa  third  parameter  L developed  by  method  8 is  valid 
only  up  to  a probability  of  suppression  evaluation  egdal  to 
1.0.  Since  for  81  mm  and  155  mm  L is  larger  than  one,  the 
function  for  the  model  should  be  decomposed.  This 
alternative  approach  to  model  suppression  is  beyond  the 
scope  of  this  paper  and  will  not  be  considered  further  at 
this  point,  although  it  is  certainly  a topic  worthy  of 
further  research. 

As  a possible  further  step,  confidence  intervals  for 
K and  £ could  be  developed. 


2»  Suppression  Time 

The  suppression  time  data,  which  were  collected  in 
seconds,  represent  the  time  for  which  an  individual  remains 
suppressed  as  a reaction  of  a single  round.  During  this 
time,  the  individual  was  unable  to  carry  out  his  mission;  he 
is  in  hiding  in  an  effort  to  survive.  According  to  the 
setup  of  the  experiment,  it  was  possible  for  the  suppressae 
to  react  to  detonations,  which  he  was  able  tc  observe  and  to 
hear  or  to  hear  only  (visual/auditory  and  auditory 
perceptors) . This  is  consistent  with  the  model  design  in 
chapter  IX. 

This  analysis  is  an  attempt  to  compare  the 
suppression  interval  gained  from  the  experiment  with  a 
certain  distribution  whose  parameters  are  estimates  of  the 
data.  The  GANNA  distribution  was  selected,  because  the  range 
cf  the  random  variable  X,  representing  the  time,  is  limited 
below  by  zero  and  by  ♦ ^ above. 

The  data  of  suppression  interval  for  the  calibers  81 
BB,  105  ma,  and  1SS  bb  are  plotted  in  histcgraas  in  appendix 
J. 
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The  estimates  A and  r for  the  Gamma  distribution 
are  derived  from  the  expectation  and  the  variance  o£  the 
data. 


E[x] 

Var(x] 


n 


1 

— J X. 

”iii  ^ 

s I 

" i=i 


A A 


r “ X • E[x]  = 


n 2 

< I ^i) 

i°l 

n , 

n I (x  -x) 


i=l 


(3.27) 


iiith  these  parameters,  the  Gamma  density  can  he  computed  in 
the  considered  range,  by  the  formula: 


f (x/X,r) 


A 


0 < X < 


0 


otherwise 


(3.28) 


?or  comparison,  the  related  density  was  computed  by 


f(x) 


frequency  In  Interval 
N ♦ interval 


(3.29) 


within  the  same  range,  where  the  frequency  in  a particular 
interval  and  the  interval  itself  is  taken  from  the 
nistograms  of  appendix  J,  and  the  constant  N is  the  sample 
size  of  the  measured  suppression  intervals. 

A.  A 

The  numerical  values  for  the  estimates  A and  r and 
the  values  for  the  sample  size  N,  the  interval  dx,  and  the 
range  are  displayed  in  the  fcllowing  table. 


Caliber 

A 

X 

A 

r 

N 

INTERVAL 

Ax 

RANGE 

SI  mm 

0.686 

2.174 

445 

0.5 

0 

to  16.15 

105  mm 

0.468 

1.374 

348 

0.4 

0 

to  13.7 

155  mm 

0.587 

2.591 

101 

0.5 

0 

to  16.75 

81 


Ibe  following  Fig  14  cospares  the  Ganna  density  (3.28) 
with  the  height  of  the  related  estinated  probability  of  the 
data  (3.  29)  . The  sketches  show  that  the  Ganna  density  for 
each  anmunition  type  respectively  underestioates  the  related 
frequency  of  the  data  because  of  its  long  tails.  In  fact, 
any  distribution  which  neets  the  requirenents  above  could  be 
taken  for  a comparison,  although  the  Ganna  distribution 
seans  to  ce  a particularly  good  choice  for  the  81  nn  caliber 
anmunition. 


Iq  cider  to  intensify  the  coaparison,  as  a next  step 
the  skewness  and  the  kurtosis  oi  the  data  and  of  the  Ganna 
distribution  with  the  estimates  above  will  be  contrasted. 

Ihese  characteristics  express  the  syooetry  of  the 
distribution  about  its  point  of  central  tendency  and  the 
relative  concentration  of  cases  at  the  center  and  along  the 
tails  of  the  distribution. 

For  the  derivation  of  the  skewness  and  the  kurtosis 
for  a GaoiBa  distribution,  we  need  the  first  four  nonents 
about  the  origin: 


« E(x]  - 


2,  r 2 -Ax  ^ ^ 

E[x  ] » / X e — X dx 

° (3.30) 


By  changing  the  variable  of  integration 
z =>  Xx 
dz  = X dx 

2 - i 


wa  find 
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^ <•  X • T dz 

r(r)  X 


(3.31) 


(3.32) 


r(r  + 2)  f 
X^r(r)  0 


T+2-1  -z 


r(r  + 2) 


dz 


tut  the  above  integral  is  equal  to  one,  hence 


^2  “ “2 


(3.33) 


K 

I 


In  a similar  fashion,  the  third  and  the  fourth  moment  can  be 
derived. 


Wd  “ (r  + 2)(r  + l)r 


(r  + 3)  (r  + 2)(r  + l)r 


y3  - M3  - 3M2  E(x]  + 2(Elx])*‘ 

M.  - M.  - AM3  E[xl  + 6u2(E[xl)^-3(E(xl)' 


(3.34) 


Converting  these  moments  to  moments  about  the  mean  by 

using  binomial  expansion; 


(3.33) 
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ue  can  coapute  the  skevness  and  the  kuctosls 


"3 

‘3  ““T 
0 


»4  4 

O 


(3.36) 


Foe  a Gaaoa  distribution  ve  get  the  following  foraulas, 
(r+2)(r-H)r  , (r+l)r  . r . , ^ 

' /F 

and  by  siaiiar  operation: 


• a. 


3r  + 6 


'4  4 

Osing  the  above  estimates  for  A and  r we  geti 


(3.37) 


a 

« ; 


<3.38) 


Ihe  formulas  for  the  skewness  and  kurtosis  derived  from  the 

data  itself  are: 

, N 

jf  I (x.  - x)- 
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(3.39) 


1 N „ 

q I 
“ 1-1  ^ 


-v2.2 


The  resulting  numerical  values  for  the  skewness  and  kurtosis 
derived  for  the  Gamma  distribution  and  from  the  data  are 
displayed  in  the  following  table: 
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81  IM 

105  Bin 

155  Bin 

Skewnesa  a, 
GAMMA 

1,35« 

1.423 

1.242 

Skewness  a. 
DATA 

2S.1S6 

1.032 

2.103 

Xurtoeis  a. 
GAMMA  ' 

5.7S9 

6.039 

5.315 

Kurtosis  a. 
DATA  * 

8.004 

1.309 

6.750 

In  discussing  the  ounerical  features  we  can  see  that  the  ' 
skewness  of  the  Gamoa' distribution  and  of  the  data  itself 
shows  an  asynmetric  right-'skewness  (positive  values). 
Although  the  values  differ  considerably  . So  is  e.  g.  the 
skewness  of  the  data  for  155  no  aloost  twice  as  big  as  the 
skewness  based  on  the  Gamaa  distribution. 

The  saoe  can  be  said  for  the  kurtosis.  There  is 
significant  difference  between  their  nuserical  values. 

This  analysis  supports  the  fact  that  the  Ganna 
distribution  can  only  be  a rough  fit  to  the  data  given,  and 
as  already  earlier  expressed,  a fit  of  another  distribution 
probably  would  have  been  as  successful  as  this  one. 


C.  VEBIfICailOM  OF  THE  MODEL 


In  ordar  to  supply  general  features  as  input  factors  for 
the  decision  process  for  military  leadership,  the  purpose  of 
this  paragraph  is  to  compute  numerical  val:.'es  for  the 
probabilities  and  rates  developed  in  chapter  TT.B.I.  and  2. 

For  the  computation  of  numerical  values,  it  is  necessary 
to  make  the  following  reasonable  assumptions: 

It  is  assumed  that  an  artillery  unit  consist  of  six 
weapons  all  either  105  mm  or  155  mm,  or  that. three 
81  ms  launchers  are  combined  to  a mission  unit. 
Hence  the  fire  rate  for  such  units  will  be 
concluded  to  be  [Uiener,p.1d9,211,213]: 

81  as  unit  with  SO,  55,  and  60  rounds/min 

105  sa  unit  with  28,  32,  and  36  rounds/min 
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I 

1 


i 


155  na  uait  with  18,21,  and  24  counds/oin 


The  standard  deviation  o for  the  density  for 
artillery  hits  will  be  considered  to  be  [FM  6-161-1, 
p.53]; 

25, 30, and  35  o for  81  Ka  unit 

/ 

20, 25, and  30  o for  105  on  unit 
30, 35, and  40  o for  155  on  unit* 

The  numbers  are  taken  from  the  field  nanual  6-161-1 
and  out  of  working  papers  of  COEC.  They  are  i 'unded 
for  convenience. 

The  parameter  H of  the  probability  of  kill  is  chcsen 
completely  arbitrarily  with  100  for  81  mm,  150  for 
105  mm,  and  200  for  155  mm  and  has  nothing  to  do 
with  experimentally  observed  values  for  the  weapons 
here  in  question. 

The  parameter  estimates  K and  c axe  taken  from 
appendix  H (tableau)  with  the  values: 


81  mm 

1 

105  mm 

155  min 

A 

k » 

1360.498 

1991.151 

2349.845 

A 

e a 

0.871 

1 

1 

0.769 

0.359 

The  followinc  figure  displays  the  conditional  probability  of 
suppression  with  the  above  estimators  for  the  parameters. 


P(S/0,r)  • exp 


r (1  - e cos  0) 


for  daO  i.e.  along  the  main  direction  of  sight. 
P(S/0  ■ 0*,r)  « exp  1^-  r >^^(1  - 


(3.40) 


(3.41) 
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Sang*  of  T t 
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R - 1991.151 
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155  ram 

K « 2349.845 
Eps  • 0.859 
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I • ♦ 
• I 


PlgUM  15  - FUNCTION  P(S/0  - 0»,r)  - cxpt-  ^ r (1  - c)J 
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For  the  above  selected  values  o£  the  standard  deviation  a, 
the  fire  rate  , and  the  paraneters  K,  e # and  H the 
values  c£  the  probabilities 


P(K)  - — — 

o /H  + 1 

and  rates 

X (with  no  killing)  » P(S)  = 

S X 

X (with  killing)  » X^  Pi,(S) 

X =■  X-  P(K) 

K ^ 


(3.43) 


(3.44) 


are  displayed  on  the  succeeding  two  figures  number  16  and 
17,  Numerical  values  for  the  rate  of  rise  are  also 

computed  and  displayed  in  figure  17.  The  collected  data 
(suppression  intervals)  were  not  distinguished  among 
individuals  as  in  the  course  of  the  foregoing  analysis, 
decause  of  this  fact  the  numerical  values  for  A are  based 
cn  the  expectation  of  duration  time  with  i ^gaid  to 
population. 


E(T] 


1 

AjP(S) 


«•  XrF(S)C 

/ e fp(L)dt  - 1 

0 ^ 


(3.43) 


It  eight  be  worthwhile  to  evaluate  the  rate  of  rise 
by  applying  the  Gamma  distributions  for  f^(t)  which  were 
discussed  in  chapter  and  contrast  them  to  the 

values  of  gained  by  the  data  Itself. 
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Hence  solving  ecuation  (3.4S)  by  having  set 


r-1 

t,(t)  - (Xt)-'  X 


(3.46) 


ue  find: 


e(t] 
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XfP(S) 


« XfP(S)t  -Xt 

/a*  ^ at) 

0 


r(r) 


(3.47) 


if  X > X^P(S)  and  by  changing  the  variable  cf  integration: 


(X  - XjP(S))t  » z (X  - X^P(S))dt  - dz 


(3.48) 


ue  get: 
E[t] 


X^P(S)[X  - X^P(S)]‘^  0 


“*  -z  r-1 
f 6 2 . 1 

J w \ dz  — 


XfP(S) 


- 1 


(3.49) 


E(t] 


X^P(S)(X  - X^P(S)l’^ 


By  using  the  earlier  developed  estimators,  we  receive  the 
mathematical  expression  for  the  rata  of  rise  based  on  a 
Gamma  di.‘’tributioQ  for  the  suppression  interval  data# 


X 

u 


X^P(S)i[X  - X^p(s)l' 


1 

XfP(S) 


(3.50) 


aemarJt:  Ihis  derivation  is  only  true  for  A > ^P(5) .which  means 
that  for  applying  this  formula,  the  fire  cate  may  not 
exceed  the  value: 


X < — n 

- P(s) 


(3.51) 


Otherwise  the  integral  (3*47)  explodes  towards  infinite 
and  the  resulting  cate  of  rise  would  ne 


(3.52) 
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which  Is  consistent  within  the  set  of  assumptions. 

Instead  of  using  the  Gamna  distribution  the  data 
themselves  were  applied  in  a second  step  to  calculate  an 
estimate  fox  and  based  on  the  number  of 
observations  o< 


(3.53) 


Ihe  approximate  value  for  the  rate  of  rise  is  then 


(3.54) 


If  killing  as  an  additional  event  is  considered,  the  rate  of 
rise  is  computed  with  the  same  formulas  (3.49)  and  (3.53) 
except  that  this  time  ^(S)  is  used  instead  of  P(S).  The 
computation  of  the  values  for  the  formulas  (3.50  ) and  (’3.54) 
was  performed  by  AFL  and  FORTRAN  respectively.  The  programs 
are  displayed  on  appendix  K. 
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Figure  17  - BATES  FOB  THB  MCCEL 
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Bow  (4)  and  (6)  and  cov  <S)  and  (7)  of  figure  17  display 
coopacdlils  values  foe  the  rate  of  rise.  As  it  was  stated 
earlier,  the  Gaama  distribution  underestimates  the  related 
freguency  of  the  data  and  hence  the  comparable  values  for 
the  rate  of  rise  in  row  (4)  are  larger  than  the  values  in 
row  (6).  The  sane  is  true  for  row  (5)  and  |7}.  They  differ 
in  general  by  10  ^ 12X.  If  we  are  willing  to  live  with  this 
fact,  the  rates  of  rise  evaluated  by  the  fitted  Gamma 
distribution  are  a good  approximation  for  the  values 
computed  by  the  data. 

In  order  to  verify  the  four  Lanchester  eguations  (2.76), 
(2.77),  (2.9fi),  and  (2.99),  presented  in  chapter  II. C.,  it 
is  necessary  to  specify  particular  rates  given  in  the  figure 
before.  In  case  of  a certain  known  composition  of  fire 
units,  specific  rates  could  be  developed  as  inputs  fer  the 
model  eguatiens. 


I- 


ay  this,  the  model  equations  receive  their  specific  shape 
and  scale,  their  general  behavior  remains  the  same,  as  can 
be  seen  in  figure  9 and  10. 


IV 


DISCUSSIQM 


Eased  oa  the  research  effort  of  this  paper  it  was  found 
that  the  phencoenon  of  suppression  as  defined  earlier',  is  a 
fflultidiiaensicnal  ptcbleia.  It  is  influenced  by 
psychological,  physiological,  and  environmental  variables. 
Ihere  exist  many  possibilities  to  ' model  dependencies  in 
general  form,  however  to  make  quantitative  statements  about 
suppression  the  modellers  have  to  restrict  their  efforts  to 
those  variables  which  are  observable  and  measurable.  Since 
the  main  objective  of  this  paper  has  been  to  establish 
models  which  are  able  to  express  relationships 
quantitatively,  the  main  thrust  has  been  to  formulate 
suppression  as  'uncticns  of  weapon  systems  and  their 
dispersions. 

For  the  evaluations  of  the  models,  a set  of  simplifying 
assumptions  was  necessary  in  order  to  guarantee  a 
mathematical  transparency.  Ihe  dependencies  developed  in 
this  thesis  postulate  some  satisfying  results  in  modelling 
suppression.  Ihe  models  reflect  sufficient  accuracy  of; 

Ihe  physionomy  of  the  human  being  and  its  resulting 
behavior  with  regard  to  suppression.  Suppression  is 
mainly  caused  by  visual  and  auditory  perception. 

Ihe  influence  of  the  weapon  systems  i.e.  their  size 
and  their  firing  capability  are  determining  factors 
for  the  amount  of  suppression. 

C£  course  the  detailed  results,  i.e.  the  estimation  of  model 
7araBete£S,  are  based  on  selected  weapon  systems  and 
scenarios.  Many  possibilities  exist  for  furth®'  'vk  in 
this  area  particularly  under  the  aspect  of  including 
extended  human  factors  components  in  the  construction  of 
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suppcessioa  oodels. 

The  authors  feel  the  paper  nay  provide  a contribution  to 
future  design  of  vargaoes  and  sinulations  as  well  as  veapon 
systens.  It  also  supports  a more  careful  analysis  of  the 
combat  situation* 
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APPENDIX  3 


DERIVATION  OF  CONSTANT  A 


The  ccnstaat  a can  be  determined  correspondingly  to  the 
following  derivation: 

Modifying  equation  (3.10)  in  CoXt,  The  Analysis  Of  Binary 
Data,  p.  33,  we  can  define  a transform  as: 


lu. 


Starting  with  the  original  model 
= P.(S/6,r)  = e ^ ^ 
and  using  the  log- transformation 


in  p, 


we  subtract 


- in 


0 

I - in  p. 

“i  / i 


and  choose  the  constant  a such  that 
ElXi  " ® 

Since  s-Binoeial  (u^P^,  n^p^Cl  - p^/) 

which  can  be  approximared  by 
- n^P(S/0,t)  + /iJ^  U 

where  U is  a random  variable  with 
EtU]  • 0 

E[U^]  - P^d-Pi) 

Inserting  S in  equation  (3) 
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(1) 


(2) 


(3) 


1 


[n  p.  + /n?  U + a 
an  = 

“iPi  J 

= E an  f 1 + f — 

L \ Vv^p^  Vi/> 

Approximating  this  by  laylor  series 


It  is  sufficient  to  consider  the  first  two  terms  of  the 
laylor  expansion  since  U and  a are  small  in  comparison  to  n. 


U a _ Ua_  _ 

j — I n.p.  2n.p.  3/2  2 « 2 2 

Pi  ^ ^ ^ ^ ’‘^i  Pi  ^’^iPi 


This  equation  has  to  be  zero  according  to  (3) 


Bence 


terms  cf  order 


terms  of  order 


: ElU]  - 0 


i ,eU--41.o 

”l  L'i  J 

a - J (1  - p^)  -0 


a - '^7  (1  - p ) 

X 

terms  of  higher  order  are  neglected. 
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APPENDIX  C 


OEBIVATiON  OF  VABIAMCE 


The  weighting  factor  Var{ln  f)  for  nethcd  7 in  chapter 
III  B is  computed  as  fellows: 

Expressirg  the  conditional  probability  of  suppression 
P(S/0,r)  as  a quotient  of  the  number  of  successes  S- and  the 
number  of  trials  n^  leads  to: 


where  the  expectation  cf  f is: 

Etf^]  =»  P^(S/8,r) 

""  (^) 

Var[£,n  f^]  * Var 

Since  the  random  number  S is  binomial  with  mean 
E(S)  • nP(S/0,r) 

and  variance 

Var(S)  - iU'(S/G,r)‘ (1  - P(S/0,r)) 

it  can  be  stated  as  a function  of  a random  variable  U with 
mean 

E(U)  - 0 

and  variance 
Var(U)  ■ i 
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Si  - ni'Pi(S/0,r)  + /ni?i(S/0,r)-(l  - Pi(S/0,r)J*U 


Inserting  this  in  the  eguatlcn  above: 


Varttn  fi]  = Var  S,n 


n^P^(S/0,r)  + /niPi(S/0,r)[l  - P^(S/0,r)].U 


/ / /Tl  - P.(S/0.r)r 

VarUn  f^l  = Var  4n  ^P^(a/e.r)-  \l  . u 


/ [1  - P,(S/6,r)] 

Vartan  fi]  = Var  an  Pi(S/a,r)  + an  ( 1 + '~n"p  (s/e/r)"  * ^ 


Since  P(S/Q#t)  is  constant  we  knew  that 

Var (an  P(S/0,r))  = 0 
which  leads  to: 

^ p (S/g^j.) 

Var(tu  - Var  Un  [l  * ^ 

Knowing  that  for  small  x's 
an(l  + x)  ~ X for  X « 
we  can  state 


/I  - P (S/0,r) 

Varlla  t^i  a Var  ^ T-p^Ts7o~.V)  ’ “ 
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Method  7 Set  Iteration 
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DATA 


^rsRn 

f 

I 

tPS 

«,t«s««e7a%} 

$*‘At  k 

•«tc»  **f  kt  *%  • 

OP  Pt  *i,t 


«*aoouiT«nto«  o.taonoiitt 


a 


t 

• I 

f 

I 

• •! 

aa  I 
«aaaa 


MPWOf  if  jlJ 
I — - 


o.naittK 

i«i7}at9«aa 

o,aaot»%ia*,a 

O.0iO1*U}'i%6 

9.0llT9>f|M«t 

i.o»Y«aoo«« 

o.s7%ar.TA«ioa 

a.TaaaiHt^a 

o.«iaaaiaa;i 

e.i}fai*t9««i 

0. ncaava^ap^ 

1. tbfa^airH 
o,a}7aai)«ki 
o.^Aaaaa^at^ 
o.ftiaaaitaai 


a.oataftafiaa^ 

*f.«»a}cy«s«ii 

*o.einv|Taa««^ 

o.naonaiaa' 

*4.s-iAn;offi 

e.t:)?d4$$ta 

••^aruaiAHc 

o.t^7»4^aiit 

*9.}a7iaQie%$ 

\ .ciaaaiiiitt 

*C.T;^^7^7l7^ 

*e.i>:v^acaoi 

”o«auaA>.'Mt 


‘DATA  105^M 


PttA 

*«.o«9<tcaiii«st« 

t 

ai»a.i477i7 

^••trricai 

6 

9.74ta4$7»4J 

crs 

o.Hiiaiiast 

SCAT  k 

fAMGC  or  It  '‘i  < 

UM9S  or  It  ^a 


9«ooo%4t«aaaoia  ‘o^aamaasa 


t 

♦a 


lomor  ifta) 


I 

9 I 

• j 

J«a  ••aa 

*«  t 

I 
I 
I 
I 
t 


t.9 


««a«| 

*• 


<).oa7aa77a^79 

Q.antoasta 

•o.iaaaiiiiaa 

9.«»7tl777ia 

fl.S7»7l?13l» 

O.An*i777a«76 

9,Sa9*CM37a 

9.03419711790 

o.o(iiaasa&a 

*9,9791707774 

*o.a74ea$)0da 

“o.^j^iaaiko 

•o.aai77jaiaa 

B*9717MOa«| 

•0.a5«7?>375a 

o.oMaaaaatyi 

t.aaeamaoi 

9.t774IU$60 

o.ae^aaaoaM 

9.1  >1719:1  III 

*o,a<)4«v«ai7i 

1 .1*>1)Ut1S 

*1  .ftioaiojn 

$.09177977994 

*1  .071»499f.l 

9,9917:»i77e»l 

•o.ftvt  ma 

*i»uiaaaa7a 

*9,917119171 

9, 4177771747 

•o.aaM^ata 

9,10iaK077« 

S?' 


DATA  ISSm 


Mro^rs 

rrn 

*a.9««n7itiAa9i 

f 

tait.iiataa 

A»«ppri(a} 

a 

ira 

•••laaivsiai 

scot  i 

»k*r.*  or  ti  *1  < 

or  ft  •!  ! 


i.ooaiaoaiuiaa  a*ttaa»«aii«f 


•i.t 


4f|l| 


•t«l 

I 

I 

I 

f 

I 

I 


r 

I. 


• !•««•« 


i.a 


•.ttaainaaii 
• »U79a«i714l 
Kivia^vva 

I,f7< 

.9.V«|99|1\M 

*i.atU9hoaia 

a,H717i.«0|9 

^•»uattia«ia 

ai«iiiKrPi 

.i«at«afaib« 

••a«i%Muau 


124 


APPENDIX  a 


THIS  PAGE  IS  HSSI  QUADlTY  FEACIICABIE 

EBPM  QQPY  gimaiTswign  th>  ra:>C;  - 


ANALYSIS  OF  PABANEISBS  K and  e 
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API  REGRESSION  PROGBANS 


KON-ITERATION  METHODS  | 

V RE  GRES  I 

Cn  2-h  36  2 pi  I 

[2]  ZC;l]-»-i^  1 

[3]  Z(:;2]^7  I 

[4]  zrx-^z  j 

C5.]  RE7A^{iWXZT+.^Z))  + .^ZT)  + ,^y  | 

C6] 

C7]  EPS-^K^BETA[.2'i  1 

C8]  rRAT-*-{BETA[.11^U)^(BKTAL22'>^V)  I 

[9]  RES-^^l-YIUT  I 

[10]  36  2 pi 

rii]  /iC;l>y  I 

[12]  AZi2l*RES  I 

7 . ■ I 


ITERATION 
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ITER  P 
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SIG-*-! 
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y**-(®((.7+(0.5x(l-p)))^/^))fSIG 
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U^UOiSIG 
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V*-V0*SIG 

[53 

Z<-  36  2 pi 

[63 

Z[;l]-ff^ 

[73 

zr;2  3'^7 

[83 

ZMZ 
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BKri3’^((r^(zr+.xz))+.xz2')+.xy 

[103 

/f'eC'D+B'rfrTri} 

[113 

EPS^X^BETXtr.2l 
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iriA  T-f-iBETIT  r 1 3 X (/o  ) + ( BET  IT  [23x70) 

[133 

RES^yC-YHAT 

[143 

36  2 pi 

[153 

A[;l3-^yC 

[163 

AC;2  3’*-/?E5 

7 


RQUBKs  SIQ  CHANGES  W ACGORANCE  TO  THE 
DIFFEREKT  ITEHATIVS  HSTHOOS. 


2 


APPliljDIX  J 


3)HIS  PAGE  IS  BEST  QUAI/lTY  PSAGIICABIE 
HSOU  copy  PUEIXISEED  TO  DD.C 


HISTOGRAMS  OF  SUPPRESSION  INTERVALS  105  vm 


****** 


■'  I 

1 I •• 


; * ^ 

{ • • 

♦—o 

««  # p4 
#«  « 

i T ^ 


j ♦— ts 

IJ 

I •; 


ki  ft  • » ^ » ft  e • 04 

« « ft  ft  ^ ft  • ft  • t « ft  ft- ft  « « • 

I ♦— <c. 

« 4 f « ft  # 

ft  ft^ft  ft  ftft  «4 

ft^lft 

ftftftft'ftv*^^ 

t 

ft  » # 

• «« 

*r: 

4ft)»««ft«C«ftft«  tfl 


. I I f 1 ( 

lftftftV''ftftftft««4«»«ftC)>ftftft«ft«ft*ft 

:«ftft4ftft««««rf»c«hftftCKCftft«ft 

!.*.*. « 


I 1 I •,  ft  ft  ft  f’  ft  ft  ft  • « «[ 

I I ! ! ’ r 

»ftft%>ft^-ftft-»«««ftft«*ftl  I 


J 4*^ 

ftVftftftHftVftCftftftftftCftftC. 
ft«t]cftfv<>»»ft#.«ftft«*4 

• ftftftftft 


. 4 « ft  *ii  «•  • 


I ••«ftftftft  ft  ftft»*»ft#ftft  ••«««*««»*  * 

* ftft^ftftlftftfttWftftftftftiftftftftftft* 


Oft#ftftft4ft4»#ft»««ftft«ftftftftllfttftftft««ftft«cft»  t 


ft«tliftftft«4ft  O 


130 


APP£MUIX  J 


HISTOGRAMS  OF  SUPPRliSSION  INTERVAIS  155  nm 


I i 


I i 


k a !•  t 

€ « «»  k o 

• « M ^ • 

I ♦•••© 


k«<«««4i#k« 
«*»»««  * 


»«)*< 


l«  kAtt#eh  ti 


If  ««ai44  «»«»44k<;4k 


I • • * « 1 1 

‘ • « % w t I 


XkX9X?»XkiK*P^77r 

«««««#. I * 


tfkK'kCiki.  * 
« • k k<'$*k*r 


• Ai«iikiik)iluekkkw«k  ;ti«i 


- . . . . 

i«i«kki«kkB#ki4iik 


i'k«kk«k«»k*kk»< 

*,  ' *i! ! » ! i»i i «« ft S« Sicii^* k k S 1 


1 •««  • k « > ■ k « k » • 4 4 9 k H • t « 4 «1  • O 


*****  ***  *1 


S . 

. » I ,• 


l->air4AflDCL. 


-lic*’  :*  vV'V.fKV" 


APPENDIX  K 


PBOGBAHS  FOR  COMPUTING  RATES  OF  RISE 


Program  for  formula  (3.54) 


FORTRAN  IV  G LEVEL  21  MAIN 


0001 

DIMENSION  TIM(500) ,SLA(20) , RETAU(20) 

0002 

N-348 

0003 

READ  (5,10)  (SLA(I) ,1=1,18) 

0004 

10 

FORMAT  (F10.5) 

0005 

READ(5,11) (TIM(I) ,I=1,N) 

0006 

11 

FORMAT (FIO. 5) 

0007 

DO  1 I«l,18 

0008 

SUM-O.O 

0009 

DO  2 Jal,N 

0010 

OEXP  ( (SLA  (I)  *T1M(  J) ) ) 

0011 

SOM»SUM+C 

0012 

2 

CONTINUE 

0013 

A»SUM/FLOAT(N) 

0014 

B-A-l.O 

0015 

ETAU*B/SLA{I) 

0016 

RETAU(I)=1.0/ETAU 

0017 

1 

CONTINUE 

0018 

* 

WRITE(6,12) (RBxAD(I) ,1=1.18) 

0019 

12 

FORMAT(1X,3F10.5) 

0020 

STOP 

0021 

END 

Program  for  formula  (3.50) 
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